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Motivated by grand unified theories and string theories we analyze the general structure of the 
neutralino sector in the USSM, an extension of the Minimal Supersymmetric Standard Model that 
involves a broken extra U(l) gauge symmetry. This supersymmetric U(l)-cxtendcd model includes 
an Abelian gauge superfield and a Higgs singlet superfield in addition to the standard gauge and 
Higgs superfields of the MSSM. The interactions between the MSSM fields and the new fields are 
in general weak and the mixing is small, so that the coupling of the two subsystems can be treated 
perturbatively. As a result, the mass spectrum and mixing matrix in the neutralino sector can be 
analyzed analytically and the structure of this 6-state system is under good theoretical control. We 
describe the decay modes of the new states and the impact of this extension on decays of the original 
MSSM neutralinos, including radiative transitions in cross-over zones. Production channels in cascade 
decays at the LHC and pair production at e + e~ colliders are also discussed. 



1 Introduction 



Adding an extra U(l)x broken gauge symmetry to the gauge symmetries of the Standard Model is well 
motivated by grand unified theories [1]. The corresponding supersymmetric extension that generalizes 
the minimal supersymmetric Standard Model (MSSM) often appears as the low energy effective theory 
of superstring theories [2]. This \J(l)x extended supersymmetric gauge theory shall henceforth be 
denoted as the USSM. 

The Higgs sector associated with the broken U(l)x gauge symmetry provides an elegant solution to 
the ll problem in supersymmetric theories [3,4]. An effective ll parameter is generated by the vacuum 
expectation value of the new singlet Higgs field S, which breaks the U(l)x gauge symmetry. This is 
the same mechanism employed by the next-to- minimal supersymmetric standard model (NMSSM) [5], 
However, the USSM possesses an additional advantage by avoiding the extra discrete symmetries of the 
NMSSM that, in the canonical version, result in the existence of domain-walls that are incompatible 
with the observed energy density of the universe. Moreover, the upper bound on the mass of the lightest 
Higgs boson of the MSSM is relaxed in the USSM due to contributions from the new singlet-doublet 
Higgs interactions and the \J(l)x D-terms [6]. Various scenarios of this type have been discussed in the 
literature, see e.g. Refs. [7,8], in which the U(l)x gauge symmetry is embedded in the grand unification 
group Eg (or one of its rank- five subgroups). 

Including the extra symmetry, the gauge group is extended to G = SU(3)c* x SU(2)^ x U(l)y x U(l)x 
with the couplings gs,g2,9Y,9x, respectively. The matter particle content in the supersymmetric theory 
includes, potentially among others, the left-handed chiral superfields Lj, Ef; Qj, Uf, Df, where the 
subscript i denotes the generation index, and the Higgs superfields H^, H u , S. The usual MSSM Yukawa 
terms Wy of the MSSM superpotential {i.e. without the ll term) are augmented by an additional term 
that couples the iso-singlet to the two iso-doublet Higgs fields: 

W = W Y + XS (H u H d ) . (1.1) 

The coupling A is dimensionless. Gauge invariance of the superpotential W under U(l)x requires the 
U(l)x charges to satisfy Qii d + QH U + Qs = and corresponding relations between the U(l)x charges of 
Higgs and matter fields. [In the following, we use Q\ = Qn d and Q2 = Qh u for notational convenience.] 
The effective ll parameter is generated by the vacuum expectation value (S) of the scalar S-field. 

Compared with the MSSM, the USSM Higgs sector is extended by a single scalar state. The neu- 
tralino sector includes an additional pair of higgsino and gaugino states, while the chargino sector 
remains unaltered. The complexity of phenomena increases dramatically by this extension but the 
structure remains transparent if the original and the new degrees of freedom are coupled weakly as 
naturally demanded [see below]. 

The supersymmetric particle spectrum of the USSM has received limited attention so far in the 
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literature [3,8-11]. In this report we attempt a systematic analytical analysis of the neutralino system, 
based on the well-motivated assumption of weak coupling between the original MSSM and the new 
additional gaugino/higgsino subsystem. In contrast to the MSSM where exact solutions of the mass 
spectrum and mixing parameters can be constructed mathematically in closed form (see e.g. Ref. [12]), 
this is not possible anymore for the supersymmetric U(l)x model in which the eigenvalue equation 
for the masses is a 6th order polynomial equation. However, analogously to the NMSSM [13], if the 
mass scales of the supersymmetric particles are set by higgsino and gaugino parameters of the order 
the supersymmetry (SUSY)-breaking scale, Msusy ~ 0(1O 3 GeV), while the interaction between the 
new singlet and the MSSM fields is of the order of the electroweak scale, v ~ O(10 2 GeV), then the 
perturbative expansion of the solution in f /Msusy provides an excellent approximation to the mass 
spectrum and yields a good understanding of the main features of the mixing matrix. 

Once the masses and mixings are determined, the couplings of the neutralinos to the electroweak 
gauge bosons and to the scalar /fermionic matter particles are fixed. Decay widths and production 
rates can subsequently be predicted for squark cascades at the LHC [14] and pair production in e + e~ 
collisions at linear colliders [15]. Of particular interest are the radiative transitions between neutralinos 
in cross-over zones, where the masses of two neutralinos are nearly degenerate. 

The report is organized as follows. In Sect. [2] we first describe the general basis of the neutralino 
sector in the USSM. Subsequently, for the naturally expected weak coupling between the MSSM and 
the new subsystem, the properties of the new higgsino and gaugino are derived in Sect.0 It is shown 
to what extent the properties of the standard neutralinos are modified. The spectrum and the mixings 
are determined analytically in a weak-coupling perturbative expansion. The neutralino masses are 
determined to second-order, whereas the mixing matrix elements are determined to first-order in the 
weak coupling. The accuracy of the perturbative results will become apparent by comparing the analytic 
approximations with the numerical solutions, thereby demonstrating that a satisfactory understanding 
of the system can be achieved. As an illustration we will study the limit in which the gaugino mass 
parameters are significantly larger than the higgsino mass parameters, where both sets of parameters are 
assumed to be much larger than the electroweak scale. A general description of the neutralino couplings 
and decay widths is given in Sect.Hl including photon transitions. We also discuss production cross 
sections in e + e~ collisions and cascade decay chains of squarks at the LHC that involve neutralinos. 
Section [5] summarizes and concludes this report. Technical details of the analytical diagonalization 
procedures for the 6x6 neutralino mass matrix for non-degenerate and degenerate levels are given in 
three appendices. 
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2 The USSM Neutralino Sector 



2.1 Supersymmetric kinetic mixing 

In a theory with two U(l) gauge symmetries, the two sectors can mix, consistently with all gauge 
symmetries, through the coupling of the kinetic parts of the two gauge bosons [16]. In the basis in 
which the couplings between matter and gauge fields have the canonical minimal-interaction form, the 
pure gauge part of the Lagrangian for the U(l)yxU(l)x theory can be written 

£gaugc = Y^ v — -^X^ X^p - X^y , (2-1) 

where the parameter sinx is introduced to characterize the gauge kinetic mixing [17]. This Lagrangian 
generalizes to 

£gau g e = ^J d2e {WyW y + WxW x + 2 sin x W Y Wx] , (2.2) 

in a supersymmetric theory, where Wy and Wx are the chiral superfields associated with the two gauge 
symmetries!*] 

The gauge/gaugino part of the Lagrangian can be converted back to the canonical form by the 
following GL(2,]R) transformation of the superfields [16,17,19]: 




1 — tan x 
1/ cos x 




(2.3) 



which acts on the gauge boson and gaugino components of the chiral superfields in the same form. The 
transformation alters the U(l)yxU(l)x part of the covariant derivative to 

D„ = d„ + igyYBp + i f—gyY tan x + —Qx) ^ (2-4) 

= 3^ + ig Y YB^ + igxQ'xB'^ ■ (2.5) 

The choice of the kinetic mixing matrix in the form given by Eq. (|2.3p is motivated by the fact that the 
hypercharge sector of the Standard Model is left unaltered by this transformation, and the new effects 
are separated in the X sector (see, e.g., Ref. [20] for an alternative choice). Consequently, the effective 
XJ(l)x charge is shifted from its original value Qx to 

Q'x = ^- 9 — Ytan X . (2.6) 

cos x 9x 

Specifically, the U(l)x charge of any field is shifted by an amount proportional to their hypercharge Y 
and the mixing parameter s'mx- Thus, as a result of the kinetic mixing, new interactions among the 



1 The normalization of the superfield W = D DV follows the conventions of Ref. [18], where V is the corresponding 
vector superfield. 
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gauge bosons and matter fields are generated even for matter fields with zero \J(l)x charge originally. 

In grand unification theories the two U(l) groups are orthogonal at the unification scale but small 
mixing [16] can be induced through loop effects when the theory evolves down to the electroweak scale. 
In string theories, kinetic mixing can be induced at the tree level [19]; however, such mixing effects must 
remain small in order to guarantee the general agreement between SM analyses and precision data in a 
natural way [21]. 



2.2 The USSM neutralino mass matrix 

The Lagrangian of the neutralino system follows from the superpotential in Eq. (jl.lj) . complemented 
by the gaugino SU(2)i, U(l)y and U(l)x mass terms of the soft-supersymmetry breaking electroweak 
Lagrangian: 

^masf n ° = —=M 2 W a W a - ^M Y YY - ^M X XX - M YX YX + h.c. 

= -^M 2 W a W a - -M\BB - -M[B'B' - M K B& + h.c. , (2.7) 

where the W a (a = 1, 2, 3), Y and X are the (two-component) SU(2)^, U(l)y and U(l)jjf gaugino fields, 
and 

Mi = My , Mi = — - ^^M Y x + M Y tan 2 X , M K = - M Y tan X • (2.8) 

cos z x cos X cos X 

In parallel to the gauge kinetic mixing discussed in Sect. 2.1, the Abelian gaugino mixing mass parameter 

Myx is assumed small compared with the mass scales of the gaugino and higgsino fields. 

After breaking the electroweak and U(l)x symmetries spontaneously due to non-zero vacuum ex- 
pectation values of the iso-doublet and the iso-singlet Higgs fields, 

i . sin 3 f \ . . cos 8 ( v \ . 1 , . 

the doublet higgsino mass and the doublet higgsino-singlet higgsino mixing parameters, 

v v 
V = *-j= and fj, x = X— , (2.10) 

are generated. The USSM neutral gaugino-higgsino mass matrix can be written in the following block 
matrix fornj 2 , 

/ M A X \ 



2 Although our initial exploratory analysis is carried out at tree-level, loop corrections can easily be included following 
the procedures of Ref. [22]. 
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where .A/f 4 is the neutral gaugino-higgsino mass matrix of the MSSM, A4 2 corresponds to the new sector 
containing the singlet higgsino (singlino) and the new U(l)-gaugino B' that is orthogonal to the bino B, 
and X describes the coupling of the two sectors via the neutralino mass matrix. More explicitly, in 
a basis of two-component spinor fields £ = (B,W 3 , H®, H®, S, B') T , the full neutralino mass matrix is 
given by [10]: 



/ 



Mt 



Mi 


—mz cp sw 
m z sp s w 





M 2 
mz cp cw 
-m z sp c w 



-mz cp sw 
mz cp cw 


-A* 



mz sp sw 
-mz sp cw 

-A* 




M K 


-fJ>x sp Q[m v C/3 

-fJ-x c/3 Q' 2 m v sp 



\ 



Q' s m s 



(2.12) 



-fi\sp -fi\cp 
V M K Q[m v cp Q' 2 m v sp Q' s m s 

where the various gaugino mass parameters Mi, M2, M[ and Mk have been defined in Eqs. (|2.7p and 
(|2.8p . Notice the absence of a diagonal mass parameter of the new singlino in contrast to the NMSSM 
where the cubic self-interaction generates this singlet mass term [13]. Two additional mass mixing 
parameters, 

m v =g x v and m s = g x v s , (2.13) 
are generated after gauge symmetry breaking and the effective charges Q' 1} Q' 2 and Q' s are defined by 

Q2 lgY x „i Qs 



Qi = 1- o — tan x, 

cos x 2g x 



Q'2 



1 9Y , 

- — tanx, 



Q's 



(2.14) 



cos x ^gx cos x 

in terms of the Qi defined below Eq. (|1.1|) . As usual, tan j3 = v 2 jv\ is the ratio of the vacuum expectation 
values of the two neutral SU(2) Higgs doublet fields, sp = sin/3, cp = cos/3, and sw,cw are the sine 
and cosine of the electroweak mixing angle 9w- 

In general, the neutralino mass matrix M§ is a complex symmetric matrix. To diagonalize this 
matrix, we introduce a unitary matrix iV 6 such that 



xi = Nl e (B,W\H d ,H u ,S,B') 



(2.15) 



where the physical neutralino states are ordered by some convention. A typical choice, motivated by 
experimental analyses, is the ordering of Xk [k = 1 5 • •? 6] according to ascending mass values. As an 
intermediate step, we shall often refer to an auxiliary convention, in which the ordering of states x^m 
denoted by primed subscripts, follows the order of the original (B, W 3 ,H^, H®, S, B') basis. 

Given the neutralino mass matrix the physical neutralino masses m? , which are real non- 
negative numbers, and the neutralino mixing matrix elements Nue can be calculated. The mass term in 
the Lagrangian is given by: 



lEK^xl + h.c., 

k=l 



(2.16) 
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The transformation of the two-component fields generates the diagonalized mass matrix for the physical 
neutralino states, 

(iV 6 )* Me (iV 6 )t = diag(mf , mf , . . . , mf ) , mf > . (2.17) 

Mathematically, this transformation is the Takagi diagonalization [23-27] of a general complex sym- 
metric matrix; see Appendix A for further details. Physically, the unitary matrix iV 6 determines the 
couplings of the mass-eigenstates xt to other particles. 

If Me is complex, then CP is violated in the neutralino sector of the theory if no diagonal matrix of 
phases P exists such that P T MeP is real. If P exists, then the neutralino interaction-eigenstates can 
be rephased to produce a real neutralino mass matrix, and the neutralino sector is CP-conserving|f| If 
Me is real, then the Takagi diagonalization of Eq. (|2.17p still applies but can easily be carried out in 
two steps. First the real symmetric matrix Me can be diagonalized by an orthogonal matrix V 6 : 

V 6 M 6 (V 6 ) T = diag(mi , m 2 , . . . , m 6 ) , (2.18) 

where the eigenvalues m& are real but not necessarily positive. The Takagi diagonalization of Me, 
which yields real non-negative diagonal mass elements, can then be achieved in a second step by taking 
mf = |mfe| and defining the unitary matrix in Eq. (|2.17p by iV 6 = (P 6 !/ 6 )*, where P 6 is a diagonal 
phase matrix with elements P^ = eV 2 6^. Here, = m^/mf = ±1 is the sign of mk, which is 
also proportional to the CP-quantum number [28] of the neutralino Xk- More precisely, the relative 
CP-quantum numbers of x,k anc ^ X?> which is the physical quantity of interest, is given by 

Although the ordering of states Xk in ascending mass values is convenient, it is often useful to 
adopt an intermediate auxiliary convention. Note that the neutralino mass matrix is easily diagonalized 
in the limit of Mk = v = (i.e., before the coupling of the MSSM with the new gaugino/singlino 
block is introduced). In this limit, Me is real after rephasing the neutralino interaction-eigenstates (if 
necessary). That is, without loss of generality, we can choose M\, M[, M<i and [i to be real in this limit, 
in which case Eq. (|2.18p yields the following mass eigenvalues: my = {M\, M2, fM, — /i, my, me'}, where 
m 5',e' = 2-^1 [l ~F a/1 + {^Q's m s/Mi) 2 \ (with m$* < me')- Away from this limit, the mass-eigenstates 
Xy will be defined such that their masses are continuously connected to the masses of the corresponding 
states in the Mk = v = limit. This defines an alternative ordering of the states x% which will be 
indicated with primed subscripts. 

We shall present a set of techniques for computing analytic approximations of the physical neutralino 
masses, m v ^ and the corresponding neutralino mixing matrix elements N®, gl . As previously indicated, 

3 In this context, the neutralino sector refers to the neutralino kinetic energy and mass terms, plus terms that couple 
the neutralinos to the gauge bosons. In this restricted sector, the neutralinos would be states of definite CP quantum 
number. Of course, it is possible to introduce CP-violating interactions through the neutralino couplings to other particles, 
e.g. matter particles of the USSM. In this case, radiative corrections could transmit these effects into the neutralino mass 
matrix. 



6 



the primed subscripts denote that these quantities refer to the physical neutralino states x?/; whose 
ordering is specified above. Of course, at the end of the computation, one can convert to an ascending 
mass ordering convention by an appropriate relabeling of the states, masses and mixing matrix elements. 



3 Small Mixing Scenarios 
3.1 General analysis 

It is well known that the MSSM neutralino mass matrix M4 can be diagonalized analytically (see, e.g., 
Ref. [13]). In contrast, the diagonalization of the new USSM 6x6 neutralino mass matrix A4q cannot 
be performed analytically in closed form. However, the case of physical interest is one in which both 
the couplings of the MSSM higgsino doublets to the singlet higgsino and to the XJ(l)x gaugino, and the 
coupling of the U(l)y and ~U(l)x gaugino singlets are weak, i.e. the elements of the 4x2 submatrix X 
in Eq. (|2.1ip are small. Then, an approximate analytical solution can be found following the procedure 
given in Appendix B. 

As an initial step, the 4x4 MSSM submatrix JH4 and the new 2x2 singlino-U(l)x gaugino submatrix 
M2 are separately diagonalized: 

M± = N 4 *M±N 4 ^ = diag(r%/,m 2 ^m3',m 4 , (3.1) 
M% = N 2 *M 2 N 2] = diag(m 5 , , m&) , (3.2) 

where the my are real and non-negative. Here we use primed subscripts to indicate that the neutralino 
states are continuously connected to the corresponding states in the Mk = v = limit, as discussed at 
the end of Sect.[2j The above procedure results in a partial Takagi diagonalization of the full neutralino 
mass matrix, 

_ / iV 4 * \ / M 4 X \ ( N 4 ^ © \ _ / M° N 4 *XN 2 ^ \ 

6= \ T N 2 * ) { X T M 2 J V <D T iV 2 t J" V N 2 *X T N^ M° ) ' 

where O is a 4 x 2 matrix of zeros. The upper left and lower right blocks of A^6 are diagonal with real 
non-negative entries, but the upper right and lower left off-diagonal blocks are non-zero. 

Performing a block-diagonalization of M.§ will remove the non-zero off-diagonal blocks while leaving 
the diagonal blocks approximately diagonal up to second order, due to the weak coupling of the two 
subsystems. That is, 

Mq = N & s M 6 N% = diag(m^ h , m^, m^, m p £, mPj", mP^) , (3.4) 

where 

_ fi / i 4X 4 - ^ 

N B ~ I 2 I x diag(e-^i' , . . . , e~ 1 ^') . (3.5) 
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A detailed derivation will be presented in Appendix B. The elements of the 4x2 mixing matrix Q are 
given by: 

RefW = ^= _ 3 , Imfij/.-/ = ^= _ 13 , 3.6) 

my — rriji my + my 

with i' = l', . . . ,4' and j' = 5', 6'. After the block-diagonalization, the upper left 4x4 and the lower 
right 2x2 blocks need not be re-diagonalized up to second order in the small mixing X between the 
blocks, but the corresponding eigenvalues are shifted m^ — > m?/ to second order in the small mixing. 
The physical neutralino masses mff are given by: 

z — ' I m.-i — m.-, my + my 



ph 

m-i ~ m 



j'=5 

■i 



P h. _ \- j [Re(iV 4 *A-iV 2 %,]2 [lm(^4^2 t) 
m£ ~ my - > <^ = = J - = — — J - — \ , [j =5,6, (3.f 

J *■ — ' m.M — m.-i m.;i -4- m..-i 



, ( ^ — my my + my 

The diagonal matrix of phases is chosen such that the m V u are real and non- negative, with the phases 
4>y given by: 

6 

<k*-Y, - ,-2 _ 2 , Re(iV 4 *AiV 2 t) t , J , Im(iV 4 *XjV 2 t), Y , [i' = 1', . . . ,4'] , (3.9) 
^ my (ml, - mi) 

4 _ 

- E - 1-2 _ 2 , Re(N 4 *XN 2 %i lm(N 4 *XN^)i,y , [f = 5', 6'] . (3.10) 
7^ m j'\ m y - m y) 

The (perturbative) Takagi diagonalization of the neutralino mass matrix M§ has now been achieved, 
with the (real and non- negative) neutralino masses given by Eqs. (|3.7p and (|3.8p . and the neutralino 
mixing matrix given by: 



?6 



N © 



N = N n I . (3.11) 

\ (D T iV 2 / K ' 



The validity of the perturbative expansion relies on the assumption tha 

Re(N 4 *XN 2 ^i,j 



<1, (3.12) 

for all choices of i' = 1', . . . , 4' and f = 5', 6'. That is, only degeneracies between the 4x4 block and 
the 2x2 block Ai 2 are potentially problematic. In particular, in the so-called cross-over zones in which 
the masses my — my exhibit a near degeneracy and the corresponding residue Re (N i *XN 2 ^)yji 7^ 0, 
mixing effects are enhanced and the analytical formalism in Appendix C must be applied. 



4 Since the my are non-negative, and by definition of order Msusy, the conditions | Im (N 4 *XN 2 )i'y /{rn^ +my ) <C 1 
are automatically satisfied. 
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3.2 The case of a real neutralino mass matrix 



We shall present numerical case studies under the assumption that the parameters of the neutralino mass 
matrix are real. The general analysis then simplifies, since a real symmetric mass matrix can always be 
diagonalized by a similarity transformation, VMV T , where V is real and orthogonal. Since some of the 
mass eigenvalues of a real symmetric matrix may be negative, we complete the Takagi diagonalization, 
N*A4N\ by introducing a suitable diagonal matrix of phases P and identifying the unitary neutralino 
mixing matrix by N = (PV)*, as indicated below Eq. (|2,18p . In this case, the (perturbative) neutralino 
mass matrix diagonalization can be performed using the three-step procedure of Ref. [13]: 

[1] Diagonalization of the submatrices A4± and M2 

In the first step, we diagonalize the (real symmetric) MSSM matrix .M4: 



The mass eigenvalues, which are real but need not be non-negative, are denoted by for i! = 1', ..,4'. 
The orthogonal diagonalization matrix V 4 is given explicitly in Ref. [12] for the most general choice 
of gaugino and higgsino mass parameters. Simple analytic forms for the neutralino mass and mixing 
matrix elements can be found in limits where either the gaugino parameters are much larger than the 
higgsino parameter or vice versa [29]. 

The exact analytic diagonalization of the new 2x2 submatrix M2 singlet higgsino-U(l)x gaugino 
submatrix M2 is straightforward. The matrix: 



Mf = V A Mi(y A ) T = diag(mi/,m 2 ',m3',m 4 /) . 



(3.13) 




(3.14) 



is diagonalized by an orthogonal rotation V 2 



MS 



M$ = V 2 M(V 2 ) T = diag(m 5 /, %) . 



(3.15) 



The eigenvalues my $ are given by 




(3.16) 



The orthogonal diagonalization matrix V 2 



is given by: 




(3.17) 



where the angle 9 S satisfies the relations: 




and 




(3.18) 
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with x = 2Q' s m s /M[. 

Two limits are of particular interest: 
(i) If m s S> then the masses and the mixing parameters are approximately given by 



m 5 / ~ —\Q' s \m s , rfiQ/ ~ \Q's\m s > and sin# s ~ sign(x)/v2 , (3.19) 

corresponding to maximal mixing due to the large off-diagonal entries in the mass matrix M.^. 

(ii) In the opposite limit, \M[\ 3> m s , and the mass eigenvalues and mixing angle are approximately 

given by 

my ~ -Qsm 2 s /M[ , m & ~ M{ + Qgm 2 /M[ , and sm9 s ~ Q' s m s /M[ . (3.20) 

This is a typical see-saw type mixing phenomenon. The heavy 6th state is a XJ(l)x gaugino-dominated 
state, whereas the 5th neutralino state is a singlet-higgsino dominated state. 

[2] Block- diagonalization of A4q 

We can now perform a block-diagonalization of Mq: 

~ ( M' D V A XV 2T \ ~ 
V 6 M 6 (V 6 ) T = V£[ „„ * „ _ n )vH T = dmg(m 1 ,,...,m 4 ,,m 5 ,,m & ), (3.21) 



V 2 X T V 4T Mn D 



where 



i4x4 — ^fifi^ 1 



— I ; , m , (3-22) 



and the elements of the real matrix are given by [cf. Eq. (|3.6p ]: 

= - , (3.23) 

my — fhji 

with i' = 1', .., 4' and f = 5', 6'. That is, the orthogonal matrix V 6 is conveniently split into the matrices 
V 4 and V 2 that diagonalize the 4x4 and 2x2 submatrices M4 and M2 respectively, and into the 
matrix that performs the subsequent block-diagonalization [13]: 

~ / V 4 O \ 
^ 6 ^ Vg . (3.24) 

After the block-diagonalization, the mass eigenvalues are shifted to second order in the perturbation X. 
The shifts are given by [cf. Eq. (|3.7|) and 



mv = m v + V [( y4Xy2T >'] 2 , [,' = 1',..,4'], (3.25) 

m,/ — mil 

^ \(y 4 xv 2T ) vv \ 2 r . . , n , . 

V ^ r^^ 1 - , J' = 5', 6'] . (3.26) 

z — ' mil — m^i 

i'=V 1 ■> 
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As expected, the eigenvalues fulfill the trace formula 

6' 

my = Mi + M 2 + M[ , (3.27) 

k'=V 

which is independent of the higgsino mass and the mixing parameters. 

The perturbative results obtained above are valid if \{V i XV 2T )ii j* / '{fny — rhj>)\ <C 1 for all possible 
choices of i! and f. In the regime of near-degeneracy, rrm — rhj>, the perturbation theory breaks down, 
and the analytic approach of Appendix C must be employed. Note that = ~ m j' is n °t a case of mass- 
eigenvalue degeneracy, so that the perturbative results obtained above should be reliable. This may seem 
to be in conflict with results of the previous subsection, since the latter corresponds to the degenerate 
case of my = rfijr, where we identify the positive masses my = \my\ m the notation of Sect. 13.11 
However, a more careful analysis reveals that the condition given by Eq. A3. 12f) does not apply, since in 
the case of opposite sign mass eig envalues, the residue Re (N A *XN 2 ^) Vjl = Re{iV A XV 2T ) Vj , = 0@ 

The higgsino doublet-singlet and the higgsino doublet-U(l)x gaugino mixings generate additional 
singlino and U(l)x gaugino components in the wave functions of the original MSSM neutralinos %° 
[i 1 = 1', ..,4'] of the size 

6' 

Y h' ~ Yl Q i'k' V k'j> [^ = l / ,-,4 / ;j' = 5 / ,6 / ] (3.28) 

k'=5' 

which is linear in the mixing parameter to first approximation as expected for off-diagonal elements. 
Reciprocally, the MSSM gaugino/higgsino components and the singlino and XJ(l)x gaugino components 
in the wave functions of x% an d X& are reduced to 

4' 

Vf, v *-J2 n l'ftfi' [i' = l',..,4';j' = 5',6'] 

v=v 

Vp w « V} v - \ {n T nV 2 ) jlk , If, k' = 5', 6'] (3.29) 
with Vp k , differing from V 2 k , only to second order in the mixing, as expected for diagonal elements. 
[3] Ensuring that the physical neutralino masses are non-negative 

The diagonalization of a real symmetric matrix by an orthogonal similarity transformation produces a 
diagonal matrix with real but not necessarily non-negative elements. Hence, some of the eigenvalues my 
will typically be negative. Defining the unitary matrix N G = (P e V 6 )*, where P 6 is a diagonal matrix 
whose k'k' element is 1 [i] if my is non-negative [negative] , the Takagi diagonalization of the neutralino 
mass matrix is achieved with non-negative neutralino masses. In particular, the unitary neutralino 
mixing matrix iV 6 * appears (instead of the real orthogonal matrix V e ) in the corresponding Feynman 
rules involving the neutralino mass-eigenstates. 

5 As in step [3] below, we identify N M = (p M V AI )* for M = 2 and 4, respectively, and {P A )~,],{P 2 )^), = -i for opposite 
sign mass eigenvalues and rhji . 
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3.3 Large gaugino mass parameters 

To illustrate the previous general discussion we shall first give a detailed parametric analysis in the limit 
in which all gaugino masses are much larger than the higgsino masses, and both sets much larger than 
the electroweak and the kinetic mixing scales, i.e. Mi, M 2 , M[ 3> fJ,, v s 3> v,Mk. All neutralino mass 
matrix parameters will be taken real. 

[1] Starting again with the diagonalization of the MSSM submatrix M4, the diagonalization matrix V 4 
defined in Eq. (|3.24p can be parameterized up to second order according to standard MSSM procedure 
(see, e.g., Ref. [12]), as 

(3.30) 

The effect of the 2x2 rotation R^u = (1 — ir 2 ) / \/2 [where t = (n , t 2 , T3) are the 2x2 Pauli matrices] 
is to shift the {34} off-diagonal elements [— fx, —fj] onto the diagonal axis [/i, —fj]. The matrix, V x , 

( — c + sw mz/Mi — c- swrmz/Mi » 




r,. 



1 c + cwmz/M 2 C-Cwmz/M 2 j 



(3.31) 



with the abbreviations c± = (cp ± s / g)/v2, removes the mixing between the blocks of the two gaugino 
and the two higgsino states. Vq and Vh rescale the gaugino and higgsino blocks themselves: 

! / s 2 w m%/Ml 

Vg ~ 12x2 - ~ 

1 \ (? w m\lMl 

l ( (l + s 2 p)M£m 2 z /2MlMi 
V H « 12x2 - 7T I , (3.32) 

1 \ (l-s 2 ^)M( / |m|/2M 1 2 M| 

with M{' 2 2 = M 2 Cyy + M^Syy. Vq and Vh relate to a diagonal form of the gaugino- higgsino mass 
matrix for large Mi 2 and fi. Their off-diagonal matrix elements are of second order and can be omitted 
consistently as they would only affect the eigenvalues at fourth order. 

The 2x2 diagonalization matrix defined in Eq. (|3.24p can be parameterized up to second order as 

V 2 „ ( 1 " Q'lm 2 s /2M? -Q' s m s /M[ \ 
~ ^ Q' s m s /M[ 1 - Ql m 2 s /2M> 2 J ' 

The 2x2 matrix V 2 generates a diagonal form of the singlino-U(l)x gaugino mass matrix for M[ 3> m s . 

After these steps are performed, the 4x4 and 2x2 mass submatrices are diagonal and the complete 
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symmetric mass matrix A4q takes the intermediate form 

/ my 




my 



my 



rri4> 



+/UAC_ -fi\c + 

\M K Q'_m v Q' + m v 








M K \ 


Q'_m v 



-H\c + Q' + m v 



m 5 / 



(3.34) 



m 6 ' / 



where, in obvious notation, zero elements of the diagonal blocks are suppressed for easier reading, and 
Q± = (Q'iC/3 ± Q' 2 sp)/V2. The diagonal elements rhk' are given by 



my 



m 7 9 



my 



111 



Z „2 



m 2 > = M 2 + ^-c w , m 4 > 
where c± is defined below Eq. (|3.3ip and 



m|Mi2 2 

u cl 

P MiM 2 + 

m|Mi2 2 

-// — — q 

^ MiM 2 ~ 



m 5 / 
m & 



M[-n K , 



M12 = M 1 c 2 v + M 2 s 2 w , 



(3.35) 



(3.36) 



The parameter can be identified with the NMSSM-type singlino mass parameter [13]. Note that 
my = fi K is small compared to all the other neutralino masses in the limit of large gaugino mass 
parameters considered in this subsection. 

[2] The block-diagonalization of the 6-dimensional intermediate matrix [Eq. (|3.34p ] can be performed by 
choosing the proper form of $7 in V®. In the limit of large gaugino mass parameters and small singlino 
mass [i K <C [i <C Mi , M 2 , M[ , the 4x2 mixing matrix Q is reduced to the simple expression 

/ M K /(Mx - M[) \ 

{i\c^/n -Q'_m v /M[ 

\ -Q' + m v /M[ J 

As a result of the block diagonalization of Eq. (I3.34|) . the mass eigenvalues are shifted according to 
Eq. (|3.25p . The resulting mass eigenvalues to the desired order are given by: 



(3.37) 



m 



my 



m 2 i 



my 



Mi + '^s 2 w + 



Mi 



Mi -M{ 



m 



Z „2 



M * + W 2 ^ 

m? z M 12 2 nlc 2 _ Q'fml 



MiM 2 



M[ 



my 



my 



m & 



m\M 12 2 
MiM 2 C " 



M|4 Q+ ml 



M[ 



Ma 

H S 2 p 

M[-fi K + 



mj{Q'l + Q'l) 
M[ 



Ml 



Mi - M[ 



(3.38) 
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Note that the sum rule given by Eq. (13.271) is satisfied. 

As expected, while the large SU(2) gaugino mass my is not affected by the singlino and the U(l)x 
gaugino, the MSSM U(l) mass my is affected by the U(l) kinetic mixing. All the higgsino states are 
modified by the interactions between the MSSM and the new subsystem. The value of my is raised by 
the interaction with the MSSM higgsinos, but remains small nevertheless. 

The mixing in the wave-functions is described by the components of O alone since the 4x4 matrix 
V 4 and the 2x2 matrix V 2 deviate from unity only to second order in the small parameters of the 
order of the SUSY scales [i = 1', ..,4']: 

2 



T/6 _ /fA m n n V 6 _ , _ Q's m s V\ 

v i'5' ~ l u ' u ' c -' c +)i' ' ~ 1 ~ ^ ,.o ~ r 

A* 

V^ i ,«-^(0,O l c^^). / , Vg, 

6 / Afg QLm v Q' + m v \ 6 Q^m^ 

6 ' ~ \Mi — M{ ' M{ ' M( ^, ' 6 ' 5 '~ M{ 



6 _ / —Mk Q\m v cp Q' 2 m vS/3 \ 6 Qgm| _ mg(Qy 2 + Q'_ 2 ) _ Afe 

yi ' \Mi-Mj' ' AfJ ' M { J., ' 6 ' 6 '~ 2Mf 2M{ 2 2(Mi- 



2(Mi - M{) 2 ' 
(3.39) 



The non-trivial mixing between two U(l) gaugino states, elements {1'6'j and {6'1'j, is generated by 
the non-zero Abelian gauge kinetic and mass mixing with non-zero Mr-. The analysis above fails when 
Mi « M{; this region of near degeneracy can be handled analytically using the results of Appendix C. 

In Eqs. (|3.38p and (|3.39j) . perturbative corrections up to second order have been included for the 
masses and diagonal mixing matrix elements, whereas only the first order corrections have been given for 
the off-diagonal mixing matrix elements. This follows the usual procedure of stationary perturbation 
theory in quantum mechanics, which associates second-order corrections to the eigenvalues with the 
first-order corrections to the wave function. Consequently, the zeros that appear in some of the matrix 
elements of V®, e ,, should be interpreted as approximate. For example, V 2 e , & and Vq, 2 > are expected to 
receive higher order perturbative corrections and hence be shifted away from zero. Nevertheless, the fact 
that the magnitude of these matrix elements are so suppressed will have some dramatic consequences 
for the behavior of the x^' an d Xe' masses m regions of near-degeneracy. 

[3] The final step is to identify N 6 = (P 6 V 6 )*, where P 6 is a diagonal matrix whose k'k' element is 1 (i) 
if my is non-negative (negative). The physical masses mP^ 1 are given by the absolute values of the m^ 
given above. The neutralino states can then be reordered in ascending (non-negative) mass if desired. 

The results of this subsection are easily generalized for the case of Mi, M 2 , M[, fi, v s ^> v, Mk- As 
long as the MSSM gaugino and higgsino parameters, M\p and \x remain significantly larger than the 
electroweak scale v, the couplings between the MSSM and the new fields, generated by X, remain weak 
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and the diagonalization of the mass matrix can still be performed analytically. However, instead of the 
approximate values 7715/ qi in Eqs. (|3.34[) and f)3.35|) the exact solutions f|3. 16|) must be used, and for V 2 
the general rotation matrix (|3.17p must be inserted. The approximation ceases to be valid at isolated 
points where X/ [rhy — fhji) is no longer a small perturbation, due to the degeneracy of mass eigenvalues 
fhi' ~ rhji. In these cross-over zones the analysis described in Appendix C must be applied. 



3.4 An illustrative example 

To illustrate the properties of the two new neutralinos and the impact of the coupling of the two 
subsystems on the original MSSM neutralinos, we study, numerically and analytically, the evolution 
of the neutralino mass spectrum and representative examples for the mixing of the particles from a 
very light new ~U(l)x gaugino across typical MSSM mass scales up to very high scales. Gauge kinetic 
mixing has only a small impact on the spectrum and it will therefore be neglected in the illustrative 
example. Throughout the evolution, including all intermediate regions, the coupling between the new 
states and the MSSM states remains weak, apart from regions of mass degeneracy. The evolution affects 
primarily the spectrum of the two new neutralino states. In the initial limit, M[ small, two medium- 
heavy degenerate states, 7715/6' ~ 0(v s ), are realized in the spectrum. At the end of the chain, M{ large, 
the spectrum is of a see-saw type, including one heavy and one nearly zero-mass state. 

As an illustrative example, we take M 2 = 1.5 TeV, m s = 1.2 TeV, [i = 0.3 TeV and M% = 0, and 
we assume the gaugino unification relation M\ = (5/3) tan 2 9\yM 2 ~ O.5M2. Also, for the numerical 
analysis in this paper, we set tan/3 = 5. We adopt the N- model charge assignments [8], 

Qi = ?=, Q2 = Q 5 = ^=- (3.40) 

2VT0 2x/l0 2x/W 



For definiteness we fix the gauge coupling at gx — 0.46, evolved from its F>& unification value of ^/5/3 gy 
down to the electroweak scale; however the results are not very sensitive to this assumption. We could 
also choose to fix Mx at its gaugino unification value under the assumption that all gaugino masses 
unify at the grand unification scale. This would correspond to a value of M[ ~ Mx = Mi = 750 GeV 
(neglecting kinetic mixing effects). However, to illustrate the structure of the system in various scenarios, 
we shall be slightly more general by allowing M[ to vary over a large range of values (0 < M[ < 5 TeV). 

To be specific, we choose the evolution with M[ 

from: M[ <C v <C fi < Mi, M 2 , v s 

to: v < n < M 1 ,M 2 ,v s < M[ . 

The evolution of the six (positive) neutralino masses^ and the values of two typical V e mixing elements, 

6 The eigenvalues 4' and 5' of the mass matrix [Eq. (|3.34[) ] are negative, while all the other eigenvalues are positive. 
Level crossing will therefore occur only between 2'-6' and 4'-5' when M[ is increased. The physical neutralino masses are 
given by the absolute values of the corresponding mass eigenvalues. 
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Neutralino Masses 




Mj [GeV] 

Figure 1: The evolution of the six neutralino masses when varying the U(l)x gaugino mass parameter 
M[ . The values used for the parameters are given in the text. The numbers with primes characterize 
the nature of the neutralino states connected with the ordering of the states when evolving from M[ = 0. 
Note that the 2' and 6' curves and the 4' and 5' curves, respectively, do not actually touch. This can be 
seen more clearly in Fig.\Q where these near intersection regions are expanded. The 1' and 5' curves, 
corresponding to opposite-sign mass eigenvalues, intersect for small M[ but affect each other only weakly. 

{5'4'} and {5'6'}, are shown in Figs. [T] and [2J The neutralino state mixings are exemplified by the V 6 
matrix elements {5'4'} and {5'6'} as representative for gaugino and higgsino mixings of the MSSM and 
the new states, as well as the mixing among the new gaugino and singlino states themselves. 

When the new U(l)x gaugino mass parameter M[ is varied from small to very large values, the 
pattern of neutralino masses evolves in an interesting way, as shown in Fig.[TJ For small M[ the set of 
parameters chosen in the previous paragraph, leads to a heavy SU(2) MSSM gaugino Xy- ^ * s followed 
by the two new states, mixed maximally in the XJ(l)x gaugino and singlino sector, an d X&- The 
fourth heaviest state is the U(l) MSSM gaugino Xi>- The lightest states are the two MSSM higgsinos x% 
and x%- If M{ is shifted to higher values, the mass eigenvalues in the new sector move apart, generating 
strong cross-over patterns whenever a mass from the new block comes close to one of the MSSM masses. 
This is realized at small M{ ~ my — Q'^m^/rhy ~ 0.91 TeV for the neutralino x% i n t ne new block and 
the SU(2) MSSM neutralino x^'i later between the new-block state x%> and the MSSM higgsino x\' f° r 
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Neutralino Mixing 

0.8 I i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i 




' 1000 2000 3000 4000 5000 

Mj [GeV] 

Figure 2: The evolution of two representative mixing matrix elements when the U(l)x gaugino mass 
parameter M[ is varied from small to large values. Large variations of the 5'4' parameter occur in the 
cross-over zone near M[ = 2.6 TeV. 

M[ « m 4 / - Q^m 2 s /rh A > « 2.68 TeV. For very large M[, x®, approaches the singlino state with a small 
mass value \fh&\ ~ Q'Jm^/M{, and x& the pure U(l)x gaugino state with a very large mass fhe' ~ M[. 

Outside the cross-over regions the approximate analytical mass spectra nearly coincide with the 
exact (numerically computed) solutions for the eigenvalues as demonstrated in Table [T] for three M[ 
values. 

The mixing pattern is more directly reflected in the elements of the rotation matrix V 6 , as shown 
in Fig.(2j For zero kinetic mixing, £5' and do n °t overlap with the U(l) MSSM gaugino, since 
V^g,, Vf, & « 0. Their overlap with the MSSM hig gsinos, Vy 4 ,, is small except in the cross-over zone. 
The mixing V®, e , between the new U(l)x gaugino and singlino states is reduced from maximal mixing 
— l/v2 for vanishing U(l)x gaugino mass parameter M[ to nearly zero mixing at asymptotically large 
M[. The moderate change in the 5'-4' cross-over zone is a reflection of the 5'4' variations by unitarity 
of the neutralino mixing matrix. 
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Table 1: Comparison between the exact and approximate neutralino masses m^o [in GeV] for three 
values of M[ . The values of the other parameters are defined in the text. 



m [GeV] 


M[ = 400 GeV 


Ml 
i 


= 2000 GeV 


Ml 
i 


= 4000 GeV 


Exact 


Appr. 


Am/m 


Exact 


Appr. 


Am/m 


Exact 


Appr. 


Am/m 


1 


294.0 


295.8 


0.6% 


294.1 


295.9 


0.6% 


211.6 


211.4 


-0.1% 


2 


302.7 


303.2 


0.1% 


301.0 


301.4 


0.1% 


294.2 


296.0 


0.6% 


3 


756.5 


755.6 


-0.1% 


380.3 


380.3 


0.0% 


304.7 


305.3 


0.2% 


4 


770.1 


770.1 


0.0% 


756.5 


755.6 


-0.1% 


756.5 


755.6 


-0.1% 


5 


1170.6 


1170.6 


0.0% 


1504.8 


1504.3 


0.0% 


1504.8 


1504.3 


0.0% 


6 


1504.8 


1504.3 


0.0% 


2379.0 


2379.0 


0.0% 


4213.9 


4213.9 


0.0% 



4 Neutralino Production and Decays 

Neutralino production rates in various channels and decay properties in various modes are affected by 
the mixing of the neutralino states and by the mass and kinetic mixings of the gauge bosons associated 
with the broken XJ(l)x and SU(2)xU(l)y gauge symmetries. 

The Z and Z' bosons can mix through kinetic coupling, as analyzed before, and mass mixing induced 
by the exchange of the Higgs fields, for example, charged under both U(l)'s. The resulting Z and Z' 
mixing is described by the mass-squared matrix 

M\ z , = ( m| A| ), (4.1) 
ZZ I Al m* 



k z " b Z' 



where the matrix elements are given by 



2 1 2 2 
m z = ^9z v , 



I' = gW {Qfc} + Q' 2 2 s}) + g 2 x v 2 sQi 
A| = \gzgxv 2 {Q'.c} - Qf 2 $) , (4.2) 

and where g 2 z = g\ + g\ . The eigenvalues of M"h z , and the Z and Z' mixing angle follow from 

m z u z 2 = \ ("4 + ™%> T \J (m 2 z - m|,) 2 + 4A|^ , 

tan20 ZZ / = -2A|/(m|, - m|) . (4.3) 

The phenomenological constraints typically require this mixing angle to be less than a few times 
10~ 3 [21], although values as much as ten times larger may be possible in some models with a light Z' 
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and reduced couplings [30]. 

For the neutralino production processes in e + e~ annihilation it is sufficient to consider the neutralino- 
neutralino-Zi^ vertices 

(x°iL\ z i\x®jL) = -gzZij cosOzz' - gxZ'a aw.6 Z z' , 



(x°L\ z 2\x iL) = +9zZij sm6 zz > ~ gxZ'n cos Z z> , (4.4) 



= +gz%ij sm.0 Z z> - gxZ'ij 

with i, j = 1, .., 6 and gz = 52/ciy ; L — ► R can be switched by substituting Z{j — ► —Z*- and Z[- — ► —Z'j*_. 
Explicitly, the couplings Zij and Z[ - are given in terms of the USSM neutralino mixing matrix N bj|3 

2kj = \ (N i3 N* 3 - N H N* 4 ) , 

Z'^ = Q[N i3 N* 3 + Q' 2 N i4 N* 4 + Q' s N l5 N* 5 . (4.5) 

Sfermion i/n-channel exchanges require the fermion-sfermion- neutralino vertices (with the fermion 
masses neglected): 

(xUhlfL) = -y/2 [g 2 (4N* 2 + ( e/ - l()N* ltw ) + gxQ'f L N* { 
(Xi L \fR\fR) = +V2 [g 2 e f twNa + g x Q'f R N i6 ] . (4.6) 

In Eq. (|4,6p the coupling to the higgsino component, which is proportional to the fermion mass, has 
been neglected. These would have to be included if one were to study, e.g., the neutralino interaction 
with the top quark and squark. 

For completeness, we also provide the fermion-fermion-Zi.2 vertices: 

{fh\Zi\fL} = -9z(l£ ~ e/svv ) cos9 Z z' - gxQ' Sh sm0 Z z> > 

(/l|-Z 2 |/l) = +gz{l[ ~ e f sw)sm0zz> ~ gxQ'f L cos9 Z z> ■ (4.7) 

When switching from L — > R in Eq. (|4,7p . the corresponding SU(2)xU(l) and U(l)x charges must be 
changed accordingly. I 3 = l{ L is the SU(2) isospin component (note that l£ R = 0), ej is the electric 
charge of the fermion / and Q'r L are the effective U(l)x charges of the left /right-handed fermions. 

The neutralino production and decay properties in the USSM model with the additional gaugino 
and singlino states depend crucially on their masses with respect to the MSSM neutralino masses. 
If they are much heavier than the other states, they will rarely be produced and so are practically 
unobservable. In contrast, if the singlino is lighter than the other states, a singlino-dominated state will 
be the lightest supersymmetric particle (LSP) into which the other neutralino states will decay, possibly 
through cascades. 

In the following subsections, we present a brief description of the general formalism of neutralino 
production and the subsequent cascade decays of the neutralinos. Once charges and mixing matrices 



7 For simplicity of notation, the USSM neutralino mixing matrix iV 6 will be denoted by N in this section. 
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are generalized to the present U(l)x theory, the phenomenological infrastructure for cross sections and 
decay widths can be copied from the MSSM. 



4.1 Singlino Production in e + e Annihilation 

The production processes of a neutralino pair in e + e~ annihilation^ 

e + e-^xh°j [*',J = l-6], (4.8) 

are generated by s-channel Z\ and Z2 exchanges, and t- and u-channel &l,r exchanges. The transition 
amplitudes, 

T ( e +e- _ x°iXf) = j Qa/3 [v(e+) ltM P a u(e~)] ["(j^Y^Cx?)] , (4.9) 

are built up by products of chiral neutralino currents and chiral fermion currents, coupled by bilinear 
"charges" Qll, Qlr etc. The four generalized bilinear charges correspond to independent helicity 
amplitudes, describing the neutralino production processes for polarized electrons/positrons [12]. They 
can be parameterized by the fermion and neutralino currents and the propagators of the exchanged 
(s)particles as follows: 

Qll = H — 2 Z \ FiL^uj H — 2 Z \ F ^L^2ij y~ ^iLj > 

s w c w s w c w s w 

Qlr = - 2 i F ^ z ii3 ~ c 2 J F ^L^2ij + -g- L * L j , 
b w L w b w c w b w 

QrL = H — 2 — 2~ F lR%lij H 2 — 5~ F2R%2ij H — -^^j ' 

Qfli? = 2 ^2 FlRZ*ij 2 ^2 F 2R Z 2ij ^ -R*-Rj . (4.10) 

s w c w s w c w s w 

The first two terms in each bilinear charge are generated by Z\ and Z2 exchanges and the third term 

by selectron exchange; D Zl 2 , DtL,R anci DuL,r denote the scaled s-channel Z\p propagators and the t- 

and it-channel left/right-type selectron propagators 

s s 
D Zi, 2 = 2 —■ fi and D (t,u)L,R = T—, 2 — ' ( 4 - n ) 

' J L,R 

with s = (p e - +p e +) 2 , t = (p e - — p.yo) 2 and u = (p_- — p^o) 2 denoting the Mandelstam variables for 

A. j A.J 

neutralino pair production in e + e~ collisions. The couplings Fn, t R of the gauge bosons Zj (i = 1,2) to 
a fermion pair are given by 

Fil = + (j(- efSw) c zz> + ^Q'f L szz>; Fir = -e/sfp c zz> + ^Q'f R s zz> , 

F 2 L = - (l[ - efSwj s Z Z' + ^Q'f L c zz r, F 2 r = +e f s 2 v s Z z> + ^Q'f R c Z Z' , (4.12) 



8 Recall that the numbering xl [* = 1, ... ,6] of the neutralinos [without primed subscripts] refers to ascending mass 
ordering. 
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where s ZZ i = sin6 ZZ :, c zz > = cos9 ZZ i, lf 3 = — 1/2 and ej = —1 for the electron charges. Finally, the 
matrices Zi^ij and the vectors Li and Ri are defined by (t\y = tanQy/) 

z m = +\ {N*Nj 3 - N i4 N*t) c zz , + (QiiViaiV/g + Q' 2 N i4 N* 4 + Q' s N i5 N* 5 ) s zz , , 

Z 2ij = ~\ {N l3 N* 3 - N i4 N* 4 ) s zz , + g {Q[N i3 N* 3 + Q' 2 N i4 N* A + Q' s N i5 N* 5 ) c zz , , 

L t = +I f 3 N i2 + ( e/ - l()t w N a + ^Q' fL N l6 , 

i?i = -e f t w Nn + —Q'tNiQ . (4.13) 

The e + e~ annihilation cross sections follow from the squares of the relevant couplings, 
o [e+e- -> = ^A P / S 2 j [l - (/if - ^ + A PS cos 2 9] Q 1 

■HfHfij Q2 + 2Ap 7 s 2 Q 3 cos 9 j d cos 9 , (4.14) 

where Sij = (l + Jjj)" 1 is a statistical factor which is equal to 1 for i ^ j and 1/2 for i = j; = m^>/yfs, 
9 is the polar angle of the produced neutralinos; and Aps = Aps(l, /i 2 , /x|) denotes the familiar 2-body 
phase space function. The quartic charges Q{ (i = 1, 2, 3) are given by products of bilinear charges: 

Qi = \ [\Qrr\ 2 + \Qll\ 2 + \Qrl\ 2 + \Qlr\ 2 ] , 
Q2 = \^[QrrQ rl + QllQI r \ , 

Q 3 = \ [\Qrr\ 2 + \Qll\ 2 - \Qrl\ 2 - \Qlr\ 2 ] ■ (4.15) 
The integration over the polar angle 9 can easily be performed analytically. 

The production cross sections for the three pairings of the two lightest neutralinos, {11}, {12} and 
{22}, are illustrated in Fig. [3] as functions of M[. For the parameter set defined in Sect. 13.41 the corre- 
sponding Z' mass is Mz 2 = 949 GeV and the ZZ' mixing angle is 6 ZZ ' = 3.3 x 10~ 3 ; these parameters 
are compatible with existing limits [21,30]. The center-of-mass energy of the e + e~ collider is set to 800 
GeV. Gf course, if Z 2 is in the reach of the collider, running on the Z 2 resonance would be the most 
natural way to explore all the facets of the new particle sector in an optimal way. 

For small values of M[ the cross section cr{xiX2l ^ s °f similar size as the MSSM prediction for 
the mixed higgsino pairs, X3/X4/ (cf- Fig.[T|), modified only by the Z' contribution. However, at and 
beyond the cross-over points with the new singlino type neutralino Xy, dramatic changes set in for 
pairs involving the lightest neutralino. Since the couplings of the mixed pair, X5/X3/, are suppressed to 
both the Z and Z' vector bosons, the cross section ^x^xl]} drops significantly. In contrast, the rising 
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Figure 3: The production cross sections for XiXij X1X2 an d X2X2 neutralino pairs in e + e collisions 
when varying the U(l)x gaugino mass parameter M[. The R and L selectron masses are chosen as 



701 GeV in this example. 



X51X51 coupling to Z' increases the cross section of the diagonal pair o~{xiXi} with rising M{. [The cross 
section for the diagonal pair c-fx^x!?} does not change as the MSSM higgsino character is modified only 
transiently in the 5'-4' cross-over zone.] 

The presence of the extra gauge boson Z2 with a mass of ~ 1 TeV alters the neutralino-pair pro- 
duction cross sections 0"{xiX?} an d ^{xljX;]} m the USSM significantly compared with the MSSM, as 
demonstrated in Tabled The production of light neutralino pairs, diagonal pairs in particular, are 
greatly enhanced although the light neutralino masses are nearly identical in the two models. 



Table 2: Comparison of production cross sections between the MSSM and the USSM. The value for 
M[ is set to zero in the USSM. For other values of M[ see Fig.\^ 



Cross Section [fb] 


^{x?x?} 


o-{xixl} 


T{X°2X°2} 


USSM 


6.5 


48.0 


6.1 


MSSM 


1.7 x 1(T 3 


67.1 


8.5 x 10~ 3 
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4.2 Neutralino cascade decays and sfermion decays 



If kinematically allowed, the two-body decays of neutralinos into a neutralino and an electroweak gauge 
bosons are among the dominant channels. The widths of the decays, xf ~^ Xj^k (& = 1)2), are 
given by 



167rm^o 

Xi 



(m| -m|o) 2 

3 —+mlo+m 2 - -2m 2 Zk 



m 



Xi 



+ 6m^om^q Re (Z^)} , (4.16) 



where Aps = Aps(l,m?o/^-0) r "x / m ~o), with 2^,- and -Z^,- defined in Eq. (|4,13p . 

Xj Xi k Xi 

Two examples, Xi ~~ * Xi^l f° r * = 6, 3, illustrate the evolution of the widths with M[ in Fig.HJ 
The neutralinos Xe an d Xi are identified with the MSSM SU(2) gaugino and the lighter of the MSSM 
higgsinos for small M{, and with the U(l)x gaugino and the singlino for large M[, respectively [cf. 
Fig.[T]. Even after Xs crosses to the U(l)x gaugino at the 2'-6' cross-over zone, the width increases 
due to an increasing phase space factor (due to the increasing mass difference) and the fact that Xe ^ ias 
a significant singlino component. However, once the state x± becomes singlino-dominated above the 
5'-4' cross-over zone, the width of the decay Xe ~~ * Xi^i drops dramatically as the mixing between the 
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Figure 4: The evolution of the neutralino decays Xi - > Xi^i ft = 6, 3/ when varying the U(l)x gaugino 
mass parameter M[. The notation follows the previous figures. 
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\J(l)x gaugino and the singlino state is strongly suppressed for large M{. The state X3 is the MSSM 
U(l) gaugino for small M{, the singlino-dominated state for moderate M[ and the heavier of the MSSM 
higgsinos for large M[. As the mass drops, even only slightly, the two-body decay X3 ~~ *■ Xi^i is 
kinematically forbidden for moderate M[. However, the mode is kinematically allowed again and its 
magnitude increases when the mass of the singlino-dominated X\ decreases sufficiently 

Similarly, the two-body decays of the charginos into a neutralino and the gauge boson are 
expected to be among the dominant channels if kinematically allowed. The widths of the decays, 
xf —> Xj , are given by 

(m 2 



r,0 w±i 



n 2 \ 1/2 
ff2 A PS 

I6irm~± 

Xi 



\w Llj \ 2 + \w Rij \ 2 



m 2 ) 



+ m 2 .± + m 2 - 



Xi 



-6m-±m^o Re(W Lij W* Rij )} , 



(4.17) 



where Aps = Xps(l,Tn'i /m 2 ± , m 2 v /m 2 ±) and the Wl,r are defined as 



Xi 1 



W Lij = U* Lil N j2 + -j=U* Li2 N j3 , W mj 



U* Ril N* 2 



1 



^-Ri2 N j4 



(4.18) 



The unitary matrices Ul and Ur diagonalize the chargino mass matrix via the singular value decomposi- 
tion [24] UrMcU^ = diag|m-±, m-± |. Explicit formulae for the chargino masses and mixing matrices 
can be found in Refs. [31,32]. 

At the LHC, sfermion decays, / — > fx® can produce complex cascades, as heavier neutralinos are 
often produced in the initial decay and subsequently decay through a number of steps before the lightest 
neutralino (which is presumably the LSP) is produced to end the chain. Thus, cascade decays are of 
great experimental interest at the LHC. The width of the sfermion 2-body decay into a fermion and a 
neutralino follows from 



r[/Wx?] 



92 A PS \ n |2 

- \9fi\ 



167rm 



m ; 



m) 



(4.19) 



where Aps = Aps(l, m 2 /m 2 ~, m 2 ^/m 2 ^), the couplings = Li and gj Ri = Ri are defined in terms of 
the neutralino mixing matrix N and the appropriate fermion charges in Eq. (|4,13l) . 



The rates for the reverse decays, neutralino decays to sfermions plus fermions, X4 
given by the corresponding partial widths^] 



//, // are 



r[x° - //] 



#2 A PS iV C 

327T m^o 

Xi 



9 



fi 



2(2: 2 

' m-o + mj 



m 



(4.20) 



As the decay rates into // and // are the same, we shall henceforth denote either of the final states by //. 
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Cascades: Sfermion and Neutralino Decays 
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Figure 5: The evolution of (a) the sfermion decays, ur — > X6 U an d ^R ~~ * X?^> an d (b) the neutralino 
decays, Xe ~^ and x® ~ > £r£> when varying the U(l)x gaugino mass parameter M[. The R-type 
slepton and R-type u-squark masses are rug = 701 Ge V and mn R = 2000 Ge V, respectively. 

with the same couplings as before, Aps = Aps(l, m 2 Jm 2 , m 2 f /rn 2 ~ ) and the color factor nL = 1,3 for 
leptons and quarks, respectively. [Analogous expressions hold for chargino decays.] 

Supersymmetric particles will be analyzed at the LHC primarily in cascade decays of some initially 
produced squark or gluino. In the U(l)x extended model, the cascade chains may be extended compared 
with the MSSM by an additional step due to the presence of two new neutralino states, for example, 

ur uxi u[Z lX i] uZ x [ti R \ uZ 1 i[t X \}. 

At each step in the decay chain, we have placed the decay products from the previous step within 
brackets. The end result of the cascade above is the final state uZ\iix\ with visible particles/jet u, 
Zi ~ Z, and two ts. 

For the parameter set introduced earlier, the partial widths involved in the cascade steps are shown 
for evolving M[ in Fig.0 The sfermion decays ur — ► XgU and £r —> x?^ are shown in the left panel and 
the neutralino decays Xg ~~ ^ Xs^i an d X5 — * ^R^ are shown in the right panel. The first step ur — ► Xgii 
corresponds to the decay of the R-type u-squark to Xg, which coincides with the MSSM SU(2) gaugino for 
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Table 3: The comparison of decay widths between the USSM and the MSSM. The state x- in the table 
denotes the second heaviest neutralino, i.e. xt in the USSM and x® in the MSSM. The value of M[ is 
set to zero in the USSM. For other values of M[ see Fig.\f^ 



Decay Width [MeV] 


T[ur -> X°iU] 




r[i R - x°A 


USSM 


130.0 


5.5 


14.1 


MSSM 


3294.6 


18.9 


15.0 



small M[ and, after the 2'-6' cross-over zone, with the U(l)x gaugino. The width increases dramatically 
before the decay is forbidden kinematically for M[ larger than 1.5 TeV. The second step Xq ~~ * x\%\ 
in this cascade chain corresponds to the decay of the MSSM SU(2) gaugino for small M[, changing to 
the U(l)x gaugino decay thereafter. The dependence of this two-body decay mode on M[ is mainly of 
kinematic nature; the decay is not allowed for M[ between ~ 0.8 TeV and ~ 1.0 TeV. Just beyond the 
2'-6' cross-over zone, it increases very sharply and keeps increasing moderately with M[ thereafter. The 
pattern for the third decay x§ ~^ ^R^ is mainly determined by the size of the U(l)y and XJ(l)x gauge 
components of the state x§- Before the 2'-6' cross-over zone the state is a U(l)x gaugino so that the 
width is large. But the width is strongly suppressed for moderate and large M[ for which x§ is a MSSM 
SU(2) gaugino with very small mixing with the two MSSM U(l) and U(l)x gauginos. The width for 
the final decay £r — > Xi^ remains moderate as the U(l) components of the Xi state are small. Beyond 
the cross-over zone, the width decreases with the suppressed U(l)x component. 

Conventional chains like q — > qx® — > q[t£\ — > q£[£xi] ma y a lso be observed in the U(l)x extended 
model. However, the partial widths in the USSM can be very different from the MSSM. As an example, 
we consider the cascade chains, in which the intermediate neutralino state Xi is the second heaviest 
neutralino, i.e. X5 m t ne USSM and X3 m t ne MSSM. As demonstrated in Tabled the width for the 
decay of ur to the second heaviest neutralino in the USSM is much smaller than in the MSSM. 

These cascade chains should only be taken as representative theoretical examples. A systematic 
phenomenological survey needs significantly more detailed analyses. 

4.3 Decays to Higgs bosons 

The USSM Higgs sector includes two Higgs doublets H u and H<i as well as the SM singlet field S [8, 
33-35]. Their interactions are determined by the gauge interactions and the superpotential in Eq. (|l,ip . 
Including soft SUSY breaking terms and radiative corrections, the resulting effective Higgs potential 
consists of four parts: 

Vh = V f + V d + U soft + AV , (4.21) 
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where the F, D and soft-breaking terms V F , Vd and V^ Q ft are given by 
V F = \X\ 2 \H U ■ H d \ 2 + A 2 |5| 2 (|i?-„| 2 + \H d \ 2 ) , 

V D = 9 f (\H d \ 2 - \H U \ 2 ) 2 + | (\H u \ 2 \H d \ 2 - \H U ■ H d \ 2 ) + 9 f (Q[\H d \ 2 + Q' 2 \H U \ 2 + Q' 5 |S| 2 ) 2 , 
F soft = m 2 \H d \ 2 + m 2 |i? u | 2 + m||S| 2 + (XA X S H u ■ H d + h.c) , (4.22) 

with H u ■ H d = H+H~ - H^H®. The structure of the F term V F is the same as in the NMSSM 
without the self-interaction of the singlet field. However the D term Vo contains a new ingredient: the 
terms proportional to g\ are D-term contributions due to the extra U(l) v which are not present in the 
MSSM or NMSSM. The soft SUSY breaking terms are collected in Vg ft- The tree- level Higgs potential 
is CP-conserving [34]. That is, one can rephase the Higgs fields to absorb the phases of the potentially 
complex coefficient XA\. Thus, without loss of generality, we will assume that these parameters are 
real. 

The term AV in Eq. (|4.2ip represents the radiative corrections to the Higgs effective potential [36]. 
The dominant contributions at one-loop are generated by top quark and scalar top quark (stop) loops 
due to the large Yukawa couplings; these terms are the same as in the MSSM. All the other model- 
dependent contributions do not contribute significantly at one- loop order [33]. Therefore, we will ignore 
these subdominant model-dependent radiative corrections in the following analysis. 

The set of soft SUSY breaking parameters in the tree-level Higgs potential includes the soft masses 
mf,m| and m 2 s and the trilinear coupling ^4^- Radiative corrections are affected by many other soft 
SUSY breaking parameters that generate masses of scalar tops and their mixings: the SU(2) and U(l) 
soft SUSY breaking scalar top masses rriQ, mjj, the stop trilinear parameter At, the supersymmetric mass 
scale Msusy and, spuriously, the renormalization scale Q. To simplify the analysis of the Higgs spectrum 
it is useful to express the soft masses m 2 , m 2 , m 2 s in terms of v s , v, tan (3 and the other parameters. The 
tree- level Higgs masses and couplings depend on four variables only: A, v s , tan (3 and A\. In the numerical 
analysis, we take 1 TeV for the new parameters, mQ,mjj, At, Q, Msusy and A\. 

Decays involving Higgs bosons can be quite different for different Higgs boson mass spectra. We 
first decompose the neutral Higgs states into real and imaginary parts as follows: 

H d = —= (v cos P + h cos (3 — H sin f3 + iA sin f3 sin ip) , 
v2 

H® = —= (v sin (3 + h sin (3 + H cos (3 + iA cos (3 sin ip) , 
V 2 

S=-=(v a + N + iA cos tp) , (4.23) 
v2 

where the CP-odd mixing angle (p is determined by tan <p = 2v s /vsin 2(3 and all the Goldstone states are 
removed by adopting the unitary gauge. Subsequently the CP-even states (h, H, N) are rotated onto 
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the mass eigenstates Hi (i = 1,2,3), labeled in order of ascending mass, by applying the orthogonal 
rotation matrix H : 



(H u H2,H 3 ) k = (h,H,N) a O* k 



(4.24) 



The resulting Higgs mass spectrum consists of three CP-even scalars, one CP-odd scalar, and two 
charged Higgs bosons. 

Generally, the width of a 2-body neutralino or chargino Xi decay to a neutralino or chargino Xj an d 
a Higgs boson 4> k (-#1,2,3 or A) is given by 



n 2 \ 1/2 
ff2 A PS 

327T771 



|C^| 2 ) +4mfcmfcRe (Cj- k C[jZ)} , (4.25) 



where Aps = Aps ( 1, m~ /mj^, rn'^/m'l.) and the left/right couplings c\j k R must be specified in each 
individual case. 



(i) For the decay of a neutralino %i to a neutralino Xj an d a scalar Higgs boson Hk, Xi ~> XjHk, 



the couplings are given by, 
C ijk(Xi ->■ X°#fc) 



(N a -Niitw)(N j3 cp-N jA 8i3)-V2— (N i3 sp+N i4 cp) N j5 

92 



+2^N l6 (C&NjaCf) + Q' 2 spN j4 ) 

92 



o H 



ik 



A 



(N i2 -N il t w )(N j 3S( 3 +N j 4Cf 3 )+V2— (N i3 cp-N i4 Sf3) N j5 

92 



+2^N l6 {Q[N j3 sp - Q' 2 c p N j4 ) 

92 



U 2k 



X 

92 



V2—N i3 N j4 - 2^Q' s N i6 N j6 



92 



ChkiXi — ► X°jHk) - C^(x° -> x°#fc) • 



(4.26) 
(4.27) 



While the first term in each of the two square brackets of Eq. (14.26f) are reminiscent of the MSSM 
couplings X^X^jh and X^X^H respectively, the other terms are genuinely new in origin, arising from the 
extra interaction terms in the USSM superpotential and the extra U(l)x gauge interactions. 

The partial widths for the kinematically allowed decays X4 5 ^ Xi#i are shown in the left panel of 
Fig. [6] as a function of M{. In the areas in which x® 5 and Xi nearly coincide with the MSSM neutralinos, 
the partial widths do not depend on M[. 

(ii) Similarly, a 2-body neutralino decay to a neutralino and a CP-odd Higgs boson, Xi ~^ Xj^> 
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Figure 6: The decay widths for x§ 4 —> X1H1 (left) an d X5 6 ~~ * Xi-^i (right) for the parameter set given 
in the text. For the purposes of example, the Higgs mass parameter Ma is set to 1.25 TeV. 



follows Eq. (|4.25p with the left/right couplings given by 



r,0 



{N i 2-N il t w ){N j3 sp-N ji cp) + V2— {N iZ cp+N iA sp) N j5 

92 



+2^N i6 (Q^Njzsp + Q' 2 N jA cp) 



A 



V2-N l3 N j4 + 2^Q' s N iG N j5 
92 92 



cos ip + (i <-> j) 



(4.28) 
(4.29) 



Again, only the first term in the square brackets is similar to the MSSM XiXjA coupling. 



The widths for the kinematically allowed decays X5 6 — > XiA are shown in the right panel of Fig. [6] 
as a function of M[ in the benchmark scenario. In contrast to the scalar case, only a few decays are 
kinematically allowed since the CP-odd scalar A is heavy. 

(iii) For completeness, we describe the decays of charginos to a neutralino and charged Higgs boson 
X~ — » X { jH~ (i = 1,2; j = 1,2, .. . ,6). These follow a similar pattern, but with the last index of the 
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N* 3 U* Lil --^ (N^+N^tw) U* Ll2 



-V2^Q' 1 N* 6 U* Ll2 --c p N* 5 U* Ll2 , (4.30) 
g 2 09 



coupling removed: 

c#(xr^%°^)=-c/3 

The same left /right couplings determine the decays of neutralinos to charginos and charged Higgs boson 
Xj — > xt H~ (i = 1,2; j = 1, 2, . . . , 6). For the parameters chosen here, the large mass of the charged 
Higgs boson allows kinematically only decays of the heavier chargino \2 t° the lightest neutralino Xi 
and H ± . 



92 

V2^Q' 2 N j6 U* Ri2 -- S pN i5 U* m2 , (4.31) 

92 92 



(iv) It is also possible for Higgs bosons to decay into neutralino/chargino states, for example the 
decays Hi — ► XiXp A — > XiX^j an d H^ — ► Xixf- Clearly this is kinematically possible only for the 
heavier Higgs states. The general form of the width for these decays fa — > XjXk (fa = Hi, A, H^), is 
given by the crossing of Eq. (|4,25|) : 



(\C[ jk \ 2 + \C* k \ 2 ) - Re{Ct jk C^)} , (4.32) 
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Figure 7: The decay widths for H3 — > X1X1 2 fl e ft) an d A ~^ X1X1 2 (right) for the same parameters as 
in Fig.\^ 
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where Aps = Aps(l, m~ /m^., ra| /m^.) and Sjk = (1 + 5jk) 1 is the usual statistical factor. The 

T /Ft 

couplings Cj -j, are related to their neutralino/chargino decay counterparts in the obvious way: 

- X?X° fe ) = cg^OcS - 7$Hi) , (4-33) 
Cij R {A - x?X°) = - xU) , (4.34) 

CJ /a (F+ - xixt) = tf/ R (xj - . (4.35) 

As an example, the Higgs boson decays H$,A — ► x? Xi 2 are displayed in Fig. [71 For small M{ the 
XyX^'-^s/^ couplings in the decays H^,A — ► X1X2 are suppressed while for large M{ the Xg/xj}'-^/-^ 
couplings are no longer suppressed. The rapid changes in the 5'-4' cross-over zone are generated by 
interference effects between the Yukawa and the gauge interaction terms. Similar interference effects, 
though less significant, occur for the decays H%,A — > AiAi near the cross-over zone. 

4.4 Neutralino radiative decays 

In the cross-over zones of the neutralino mass eigenvalues, the gaps between the neutralino masses 
become very small. As a result, standard decay channels are almost shut and photon transitions between 
neutralino states [37] become enhanced. These photon transitions are particularly important in the 
cross-over zone 4'-5' at M[ ~ 2.6 TeV [cf. Fig.[TJ. The proximity of the two heavier states to the lightest 
neutralino dramatically reduces the rates of all other decay modes so that the radiative decays 

X 2,X°3^Xi+7 and X3 ^ X2 + 7 , (4-36) 

become non-negligible modes. Of course, also the 7 transitions are phase-space suppressed in cross-over 
zones but less strongly than the competing standard channels due to the vanishing photon mass, even 
for 3-particle decays into a lighter neutralino and lepton- or light-quark pair. 

The effective couplings g^o^o^ in the partial decay widths 

#?o- (m~o - m 2 ) 3 
T[X° - XH = *' 5 Xj , (4-37) 

x" 

are of magnetic or electric dipole type depending on the relative CP quantum numbers of x^ an d 
X?- The couplings are generated by triangle graphs of sfermion/fermion, chargino/VF-boson and 
chargino/charged Higgs-boson lines. The sum of all two-point graphs associated with the photon line 
and attached to the neutralino legs by a Z-boson line vanish in the non-linear .R-gauge [37]. The 7 
transition amplitudes are finally complex combinations of mixing matrix elements with reduced triangle 
functions. 

For the 7 transitions of Eq. (|4.36p . the partial widths are displayed in Fig. [8] for the set of parameters 
chosen earlier. In this example the three lightest neutralino states are predominantly of higgsino type 
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Figure 8: Comparison between photon transitions [full lines] and ordinary lepton-pair decays [dashed 
lines] of neutralinos in the 3' -4' -5' cross-over zone of the reference point. The U(l)x gaugino mass 
parameter M[ is varied between 2 and 4 TeV for comparison with the side-band wings; the charged 
Higgs mass, affecting the photon loop couplings is chosen 1.25 TeV and the trilinear parameters are set 
to 1 TeV. 

so that decays to lepton pairs are allowed through couplings mediated by virtual Z-bosons which in 
fact are the dominant modes. Therefore, for illustrative comparison, the partial widths for standard 
electron-pair decays Xi — > \j e + e~ are also shown. Evidently the branching fractions for radiative 
decays in general drastically change in the complex 4'-5' cross-over zone near M[ ~ 2.6 TeV compared 
with the side-band wings. 

5 Summary and Conclusions 

In this report we have investigated the neutralino sector in the XJ(l)x extension of the minimal Super- 
symmetric Standard Model, as suggested by many GUT and superstring models. The extended model 
has attractive features which solve several problems of the MSSM. It provides a natural solution of the 
//-problem without creating cosmological problems. The upper limit of the light Higgs mass is somewhat 
increased and the growing fine tuning in this sector is reduced. 

While the MSSM neutralino sector is already quite complex, the complexity increases dramatically in 
the extended model due to two additional degrees of freedom. However, the small coupling between the 
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original four MSSM states and the two new states offers an elegant analytical solution of this problem 
within a perturbative expansion. We have worked out this solution in detail for the mass spectrum and 
the mixing of the states. 

The expansion in the parameter v/M^vsy, the ratio of the electroweak scale v over the generic 
supersymmetry-breaking scale MsusY; leads to an excellent approximation of the exact solutions. Even 
in the cross-over zones in which two mass eigenvalues are nearly degenerate, proper adaption of the 
analytical formalism provides an accurate description of the system. Thus, in the limit in which MgusY 
is sufficiently above the electroweak scale, the neutralino system of the U(l)x-extended supersymmetric 
standard model is under good analytical control and its features are theoretically well understood. 

A few examples of mass spectra, widths for cascade decays at LHC, decays to Higgs bosons, photon 
transitions and production cross sections in e + e~ collisions illustrate the characteristic features of the 
model. 
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Appendix A: Takagi diagonalization of a complex symmetric matrix 

In quantum field theory, the most general neutral fermion mass matrix, M, is complex and symmetric. 
To identify the physical eigenstates, M must be diagonalized0 However, the equation that governs 
the identification of the physical fermion states is not the standard unitary similarity transformation. 
Instead it is a different diagonalization equation that was discovered by Takagi [23], and rediscovered 
many times since [24] § Despite this illustrious history, the mathematics of the Takagi diagonalization is 
relatively unknown among physicists. Thus, in this appendix we present a self-contained introduction to 
the Takagi diagonalization of a complex symmetric matrix. After presenting some background material 

10 An alternative method — the standard diagonalization of the hermitian matrix M* M, which is commonly advocated in 
the literature, fails to identify the physical states in the case of mass-degenerate fermions, as noted below Eq. (|A,3[) . 

11 Subsequently, it was recognized in Ref. [38] that the Takagi diagonalization was first established for nonsingular complex 
symmetric matrices by Autonne [39]. 
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and a constructive proof of Takagi's result, we provide, as a pedagogical example, the explicit Takagi 
diagonalization of an arbitrary 2x2 matrix. The latter will be particularly useful for considering cases 
in which there is a near-degeneracy in mass between two of the neutral fermions. 



A.l General analysis 

Consider a system of n two component fermion fields £ = (£1 , £2 j • • • > £,n) T , whose physical masses are 
governed by the Lagrangian 

£mass = ^ T M£ + h.C. (A.l) 

In general, the mass matrix M is an n x n complex symmetric matrix. In order to identify the physical 
masses m; and the corresponding physical fermion fields Xii one introduces a unitary matrix U such 
that £ = Ux and demands that £ T M £ = rriiXiXi- This corresponds to the Takagi diagonalization of 
a complex symmetric matrix^ which is governed by the following theorem [23,24]: 

Theorem: For any complex symmetric nxn matrix M, there exists a unitary matrix U such that 1^1 

U T M U = M D = diag(mi, m 2 , . . . , m n ) , (A.2) 

where the m k are real and non-negative. 

The uik are not the eigenvalues of M. Rather, the m k are the so-called singular values of the 
symmetric matrix M, which are defined to be the non-negative square roots of the eigenvalues of M^M. 
To compute the singular values, note that: 

U^M^MU = M 2 D = diag(mf, mj, . . . , m 2 n ) . (A.3) 

Since M^M is hermitian, it can be diagonalized by a unitary similarity transformation. Although U can 
be determined from Eq. (|A.3P in cases of non-degenerate singular values, the case of degenerate singular 
values is less straightforward. For example, if M = ^ J V the singular value 1 is doubly-degenerate, 
but Eq. (|A.3[) yields U^U = 1-2x2, which does not specify U. Below, we shall present a constructive 
method for determining U that is applicable in both the non-degenerate and the degenerate cases. 

Eq. (|A.2|) can be rewritten as MU = U*Me>, where the columns of U are orthonormal. If we denote 
the kth column of U by v^, then, 

Mv k = m k v* k , (A.4) 



12 If U — N* , we obtain the form of the Takagi diagonalization used in Eqs. (|2.17|l and (|B.2[) . 

13 In Ref. [24], Eq. (|A.2f) is called the Takagi factorization of a complex symmetric matrix. We choose to refer to this 
as Takagi diagonalization to emphasize and contrast this with the more standard diagonalization of normal matrices by 
a unitary similarity transformation. In particular, not all complex symmetric matrices are diagonalizable by a similarity 
transformation, whereas complex symmetric matrices are always Takagi-diagonalizable. 
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where the to& are the singular values and the vectors Vk are normalized to have unit norm. Following 
Ref. [25], the i>& are called the Takagi vectors of the symmetric complex n x n matrix M. The Takagi 
vectors corresponding to non-degenerate non-zero [zero] singular values are unique up to an overall 
sign [phase]. Any orthogonal [unitary] linear combination of Takagi vectors corresponding to a set of 
degenerate non-zero [zero] singular values is also a Takagi vector corresponding to the same singular 
value. Using these results, one can determine the degree of non-uniqueness of the matrix U. For 
definiteness, we fix an ordering of the diagonal elements of Mb- If the singular values of M are distinct, 
then the matrix U is uniquely determined up to multiplication by a diagonal matrix whose entries 
are either ±1. If there are degeneracies corresponding to non-zero singular values, then within the 
degenerate subspace, U is unique up to multiplication on the right by an arbitrary orthogonal matrix. 
Finally, in the subspace corresponding to zero singular values, U is unique up to multiplication on the 
right by an arbitrary unitary matrix. 

We shall establish the Takagi diagonalization of a complex symmetric matrix by formulating an 
algorithm for constructing U. A method will be provided for determining the orthonormal Takagi 
vectors Vk that make up the columns of U. This is achieved by rewriting the n x n complex matrix 
equation Eq. (jA.4j) [with m real and non-negative] as a 2n x 2n real matrix equation [40]: 




FieM 
ImM 




where m > . (A. 5) 



Since M = M T , the 2nx2n matrix M s defined by Eq. (|A.5|) is a real symmetric matrix. In particular, Ms 
is diagonalizable by a real orthogonal similarity transformation, and its eigenvalues are real. Moreover, 
if m is an eigenvalue of M s with eigenvector (Reu, Imv), then — m is an eigenvalue of M s with 
(orthogonal) eigenvector (— Imv, Rev). This observation proves that Ms has an equal number of 
positive and negative eigenvalues and an even number of zero eigenvalues l 14 l Thus, Eq. (|A.4p has been 
converted into an ordinary eigenvalue problem for a real symmetric matrix. Since m > 0, we solve the 
eigenvalue problem M s u = mu for the eigenvectors corresponding to the non-negative eigenvalues f^l It 
is straightforward to prove that the total number of linearly independent Takagi vectors is equal to n. 
Simply note that the orthogonality of (Re«i, Imv\) and (— Imv\ , Revi) with (Ret>2 , Imt^) implies 
that v\v2 = 0. 

Thus, we have derived a constructive method for obtaining the Takagi vectors v^. If there are 
degeneracies, one can always choose the Vk in the degenerate subspace to be orthonormal. The Takagi 
vectors then make up the columns of the matrix U in Eq. (|A.2p . A numerical package for performing 



14 Note that (—Imv , Rev) corresponds to replacing Vk in Eq. (|A.4|I by ivu- However, for m < these solutions are not 
relevant for Takagi diagonalization (where the nit are by definition non-negative). The case of m — is considered in 
footnote 1151 

lo For m = 0, the corresponding vectors (Reu, Imv) and (—Imv, Rev) are two linearly independent eigenvectors of 
Ms; but these yield only one independent Takagi vector v (since v and iv are linearly dependent). See footnote 1 141 
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the Takagi diagonalization of a complex symmetric matrix has recently been presented in Ref. [27] (see 
also Refs. [25,26] for previous numerical approaches to Takagi diagonalization). 



A. 2 Example: Takagi diagonalization of a 2 x 2 complex symmetric matrix 

The Takagi diagonalization of a 2 x 2 complex symmetric matrix can be performed analytically!"^! The 
result is somewhat more complicated than the standard diagonalization of a 2 x 2 hermitian matrix by 
a unitary similarity transformation. Nevertheless, the corresponding analytic formulae for the Takagi 
diagonalization will prove useful in Appendix C in the treatment of nearly degenerate states. Consider 
the complex symmetric matrix: 

Ca c \ , . , 

c 6 ) ■ (A 6> 

where c / and, without loss of generality, \a\ < \b\. We parameterize the 2x2 unitary matrix U in 
Eq.COD by [41]: 

/ cos 9 e 1 ^ sin 9 
U = VP= _ _ _ , (A.7) 



-e r sintf cos l 







where < 9 < tt/2 and < a, /3 , (ft < 2tt. However, we may restrict the angular parameter space 
further. Since the normalized Takagi vectors are unique up to an overall sign if the corresponding 
singular values are non-degenerate and non-zero0 one may restrict a and fti to the range < a , (3 < tt 
without loss of generality. Finally, we may restrict 9 to the range < 9 < 7r/4. This range corresponds 
to one of two possible orderings of the singular values in the diagonal matrix . 



Using the transformation (|A.7|) . we can rewrite the Takagi equation (|A.2p as follows: 

CTl \ 
G2 ) ' 




(A.« 



where 

cti = mi e 2ia , and a 2 = m 2 e 2l/3 , (A.9) 

with real and non-negative m^. Multiplying out the matrices in Eq. (|A.8|) yields: 

ai =a- ce-^tg = be~ 2i ^ - ce^tg 1 , (A.10) 
a 2 = b + ce i % =oe 2 ^ + ce i V 1 (A.ll) 



16 The main results of this subsection have been obtained, e.g., in Ref. [27]. Nevertheless, we provide some of the details 
here, which include minor improvements over the results previously obtained. 

17 In the case of a zero singular value or a pair of degenerate of singular values, there is more freedom in defining the 
Takagi vectors as discussed below Eq. (|A,4[) . These cases will be treated separately at the end of this subsection. 
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where tg = tan#. Using either Eq. (1A.10I) or (jA.llh . one immediately obtains a simple equation for 
tan 26 = 2(t~ 1 - tg)- 1 : 

2c 



tan 29 



(A.12) 



be- i( t> — aeW 

Since tan 26 is real, it follows that be* e~ 1 ^ — ac* e 1 ^ is real and must be equal to its complex conjugate. 
The resulting equation can be solved for e 2i< ^\ 

be* + a*c 



or equivalently 



3 2i</> 



J<t> 



b*c + ac* 
be* + a*c 



(A.13) 



(A.14) 



| be* + a*c\ 

The (positive) choice of sign in Eq. (1A.14|) follows from the fact that tan 26 > (since by assumption, 
< 6 < 7r/4), which implies < c*(be~ % ^ — ae 1 ^) = |c| 2 (|6| 2 — |a| 2 ) after inserting the results of 
Eq. (|A.14j) . Since \b\ > \a\ by assumption, the asserted inequality holds as required. 



Inserting the result for e 1 ^ back into Eq. (|A.12|) yields: 

2\bc* +a*c\ 



tan 2(9 



(A.15) 



One can compute tan^ in terms of tan 26 for < 6 < 7r/4: 



tm\6 



1 



tan 26 



y/l + tan 2 26 - 1 



|&| 2 + V(|&| 2 -|o| 2 ) 2 + 4|bc* + a*c| 2 
2\bc* + a*c\ 

2\bc* + a*c\ 
\a\ 2 + ^(l&l 2 - |a| 2 ) 2 + 4|6c* + a*c\ 2 ' 



(A.16) 



(A.17) 



Starting from Eqs. (|A. 10|) and (jA.llj) . it is now straightforward, using Eqs. (jA.14|) and (|A.16|) . to compute 
the squared magnitudes of a^'- 

l r, 



2 i i2 
m k = Wk\ 



\a\ 2 + \b\ 2 + 2\c\ 2 =f V(l&l 2 -|a| 2 ) 2 + 4|6c*+a* C i 2 



(A.18) 



with | cTi| < | o"2 1 - This ordering of the is governed by the convention that < 6 < 7r/4 (the opposite 
ordering would occur for 7r/4 < 6 < vr/2). Indeed, one can check explicitly that the \(Jk\ 2 are the 
eigenvalues of M^M, which provides the more direct way of computing the singular values. 

The final step of the computation is the determination of the angles a and j3 from Eq. (|A.9p . Inserting 
Eqs. (IAT41) and (lA~T7l) into Eqs. (lATOl) and (jATTTjl . we end up with: 



a = iarg{a(|6| 2 -|ai| 2 )-6*c 2 }, 
0= iarg{6(|cT2| 2 - |a| 2 ) + a*c 2 }. 



(A.19) 
(A.20) 



37 



If det M = ab — c 2 = (with M^O) , then there is one singular value which is equal to zero. In this 
case, it is easy to verify that a\ = and |<72| 2 = TV (M^M) = \a\ 2 + \b\ 2 + 2|c| 2 . All the results obtained 
above remain valid, except that a is undefined [since in this case, the argument of arg in Eq. (IA.19j) 
vanishes]. This corresponds to the fact that for a zero singular value, the corresponding (normalized) 
Takagi vector is only unique up to an overall arbitrary phase [cf. footnote [T7] . 

We provide one illuminating example of the above results. Consider the complex symmetric matrix: 

*=(;_;)■ (A*) 

The eigenvalues of M are degenerate and equal to zero. However, there is only one linearly independent 
eigenvector, which is proportional to (1 , i). Thus, M cannot be diagonalized by a similarity transfor- 
mation [24]. In contrast, all complex symmetric matrices are Takagi-diagonalizable. The singular values 
of M are and 2 (since these are the non-negative square roots of the eigenvalues of M*M), which are 
not degenerate. Thus, all the formulae derived above apply in this case. One quickly determines that 
9 = 7r/4, (f> = it/2, (3 = 7r/2 and a is indeterminate (so one is free to choose a = 0). The resulting 
Takagi diagonalization is U T MU = diag(0 , 2) with: 

j . (A.22) 

This example clearly indicates the distinction between the (absolute values of the) eigenvalues of M and 
its singular values. It also exhibits the fact that one cannot always perform a Takagi diagonalization by 
using the standard techniques for computing eigenvalues and eigenvectors!"^! 

We end this subsection by treating the case of degenerate (non-zero) singular values, which arises 
when be* = —a*c. Special considerations are required since not all the formulae derived above are 
applicable to this case [cf. footnote [IT]. The condition be* = —a*c implies that |o| = \b\, so that 
|cri | 2 = \a 2 \ 2 = \b\ 2 + \c\ 2 . After noting that a/c = -b*/c*, Eq. (|A"l2l) then yields: 

tan 29 = [Re(6/c) c + Im (b/c) s^]" 1 , (A.23) 

where = cos (j> and = sin eft. The reality of tan 29 imposes no constraint on (/); hence, (/) is 
indeterminate [a fact that is suggested by Eq. (|A.14[) ] . The same conclusion also follows immediately 
from Eq. (fP|) . Namely, if M D = ml 2X 2, then (UO) T M(UO) = O t M d O = M D for any real orthogonal 
matrix O. In particular, (f> simply represents the freedom to choose O [see, e.g., Eq. (|A.28|) ]. Since (f> 
is indeterminate, Eq. (|A.23P implies that 9 is indeterminate as well. In practice, it is often simplest 






18 For real symmetric matrices M, one can always find a real orthogonal V such that V T MV is diagonal. In this 
case the Takagi diagonalization is achieved by U = VP, where P is a diagonal matrix whose kk element is 1 [i] if the 
corresponding eigenvalue ra^ is positive (negative). Of course, this procedure fails for complex symmetric matrices [such 
a.s M in Eq. (|X2lj) ] that are not diagonalizable. 
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to choose a convenient value, say = 0, which would then fix 9 such that tan 26 = [He (b/c)] . For 
pedagogical reasons, we shall keep (f> as a free parameter below. 

Naively, it appears that a and (5 are also indeterminates. After all, the arguments of arg in both 
Eqs. (|A.19|) and (|A.20|) vanish in the degenerate limit. However, this is not a correct conclusion, as the 
derivation of Eqs. (|A.19jl and (IA.20[) involve a division by 1 6c* + a*c\, which vanishes in the degenerate 
limit. Thus, to determine a and (3 in the degenerate case, one must return to Eqs. (lA.lOj) and (jA.lip . 
A straightforward calculation [which uses Eq. (|A.23[) ] yields: 



(72 

C 



c* 



which implies 



a + (3 = arg c ± 



7T 



(A.24) 



(A.25) 



Note that separately, a and (5 depend on the choice of cj) (although drops out in the sum). Explicitly, 

-i<j> 



0"! 



-ce 



1 + 



Re (b/c) + Im (b/c 



+ i 



s^Re (b/c) — c^Im (b/c) 



(A.26) 



a 2 



c e 



ii)> 



1 + 



c^Re (b/c) + s^Im (6/c) 



Re (6/c) — Im (6/c) 



(A.27) 



One easily verifies that Eq. (|A.24p is satisfied. Moreover, using Eq. (|AI)|) . a and /3 are now separately 
determined. 



We illustrate the above results with the classic case of M = (® qJ. In this case M^M = l2x2 5 so 
U cannot be deduced by diagonalizing M^M. Setting a = b = and c = 1 in the above formulae, it 
follows that 9 = 7r/4, a x = -e'^ and a 2 = e^, which yields a = -(cb ± vr)/2 and = cf)/2. Thus, 
Eq. (IAT7|) yields: 

1 / 1 \ ( ±ie irt >/ 2 \ if ±ie i(t> l 2 e^l 2 \ 



U 



V2 



1 







-i(j>/2 



V2 \ ^ie-W e-^l 2 



1 

V2 



(A.28) 



t 1 \ / ±cos(>/2) sin(0/2) \ 

-i 1 J \ Tsin(0/2) cos(0/2) / ' 

which illustrates explicitly that in the degenerate case, U is unique only up to multiplication on the 
right by an arbitrary orthogonal matrix. 



Appendix B: The small-mixing approximation 

The 6x6 USSM neutralino mass matrix of Eq. (|2.12p cannot be diagonalized analytically in general. 
However, simple analytical expressions for masses and mixing parameters can be found, similarly as in 
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the NMSSM, by making use of approximations based on the natural assumption of small doublet-singlet 
higgsino, doublet higgsino-U(l)x gaugino mixing and kinetic gaugino mixing, i.e. for a large SUSY scale 
compared to the electroweak scale. 

In this appendix, we provide details of the neutralino mass matrix diagonalization in the small mixing 
approximation, in which the weak coupling between two off-diagonal matrix blocks can be perturbatively 
treated. For mathematical clarity, we present the solution for a general complex (JV + M) x (N + M) 
symmetric matrix in which the N x N and M x M submatrices are coupled weakly so that their mixing 
is small: 

/ M N X NM \ 
M N+M = (B.l) 
V X nm Mm J 

To obtain the corresponding physical neutralino masses, one must perform a Takagi diagonalization of 

( N N+My Mn+m (Ar JV+A/ )t = diag(mi , , m2 , , . . . , m N , +M ,) , my > , (B.2) 

where N N+M is a unitary matrix^ The Takagi diagonalization of a general complex symmetric matrix 
is described in Appendix A. The non-negative my are called the singular values of M, which are defined 
as the non-negative square roots of the eigenvalues of M^M. 

In Eq. (jB.ip . Mn and Mm are N x N and M x M complex symmetric submatrices with singular 
values generally of the SUSY scale, Msusy- Xnm is a rectangular N x M matrix whose matrix elements 
are generally of the electroweak scale. Assuming that the electroweak scale is significantly smaller than 
Msusy, one can treat Xnm as a perturbation as long as there are no accidental near-degeneracies 
between the singular values of Mn and Mmi respectively. (The case of such a near-degeneracy is the 
subject of Appendix C.) The diagonalization of Mn+m can be performed using the following steps. 

[1] In the first step, we separately perform a Takagi diagonalization of Mn and Mm- 

M N = N N *M N N N ^ = diag(mi/, ...m N >), (B.3) 
M D M = N M *M M N An = diag(mjv, +1 ,, . . . ,m N '+M') , (B.4) 

where the my are real and non-negative. The ordering of the diagonal elements abovj^l is chosen 
according to some convenient criterion (e.g., see the discussion at the end of Sect. [21) Analytical expres- 
sions can be obtained for the singular values and the Takagi vectors that comprise the columns of the 

19 In Eq. ()B.2|) . we use primed subscripts to indicate that the corresponding neutralino states are continuously connected 
to the states of the unperturbed block matrix, diag(A / f^ , Mm)i where the diagonal matrices Mn and Mm are defined 
in Eqs. pO]) and pi} . 

20 When N and M are used in subscripts and superscripts of matrices, they refer to the dimension of the corresponding 
square matrices. For rectangular matrices, two subscripts will be used. 
21 See footnote [191 
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corresponding unitary matrices N N and N M for values of N, M < 4 [12]. 
Step [1] results in a partial Takagi diagonalization of Mn+m '■ 



M 



N+M = 



N N * 

qT N M * 



Xnm Mm 



T N AI t 



N N *X NM N M ^ 

-T7 D 

M M 



(B.5) 



where is an N x M matrix of zeros. The upper left and lower right blocks of Mn+m are diagonal 
with real non-negative entries, but the upper right and lower left off-diagonal blocks are non-zero. 

[2] The ensuing (N + M) x [N + M) matrix, Mn+m, can be subsequently block-diagonalized by 
performing an (N + M) x (N + M) Takagi diagonalization of Mn+m- Since the elements of the off- 
diagonal blocks of M. n+m are small compared to the diagonal elements my , we may treat Xnm as a 
perturbation. More precisely, Xnm can be treated as a perturbation if: 

Re(N N *X NM N M ^;, 



< l 



(B.6) 



for all choices of i' = 1', 
computation below. 



mi' — mf 

N' and j' = N' + 1' . . . , N' + M' . This condition will be an output of our 



The perturbative block-diagonalization is accomplished by introducing an {N + M) x (N + M) 
unitary matrix: 



A* 



1 



NxN 



B 



-fit 



<-MxM 



±fitfi 



(B.7) 



where fi is an N x M complex matrix that vanishes when Xnm vanishes (and hence like Xnm is 
perturbatively small). Note that TJbJ^b = ^(n+m)x(n+m) + 0(fi 4 ) which is sufficiently close to the 
identity matrix for our purposes. Straightforward matrix multiplication then yields: 



AT 



M 



B 



B 



B T M 



D 
M 



M'j^ + 0{BVL A ) 
B T + M M 9t -Q T M N + 0{BQ 2 ) 



b + n*M M - M N n + o(Bn 2 ) 
M'jfi + o(Bn 3 ) 



(B.f 



where 



B = N N *X NM N Mjf 



M' N J =M N + 
M'^ = M M - 



n*B T + hn*M M tf - ±M^fifit + transp 



D. 



n T B 



±fi T mZsi + ±M^fitfi + tr 



ansp 



(B.9) 
(B.10) 

(B.ll) 
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and "transp" is an instruction to take the transpose of the preceding terms inside the bracket. For a 
consistent perturbative expansion, we may neglect all terms above that are hidden inside the various 
order symbols in Eq. (|B.8h . Hence, a successful block-diagonalization is achieved by demanding that 

B = M°Sl - SI*Mm . (B.12) 

Inserting this result in Eqs. (jB.lOp and ([B.lip and eliminating B, we obtain: 

n*M^ ~ Mfarf + transpj , (B.13) 

n T M^n-M^n + tTajasp\ . (B.14) 



KAlD _ -T1 D _ 1 
JVl N — JVl N 2 



M>£ = M D M -\ 



The results above simplify somewhat when we recall that and Ai^i are diagonal matrices [see 
Eq. (jR3j) and (|B~4|) ]. Taking the real and imaginary parts of the matrix elements of Eq. (|B.12|) yields 
two equations for the real and imaginary parts of n^-: 

Refi,, v = = ^— , Imfi^v = = ^— , B.15 

my — rriji my + my 

with i' = 1', . . . ,N' and f = N' + 1' . . . , N' + M' . Since the Qyji are the small parameters of the 
perturbation expansion, it follows that | Re By j> /(my — rnj>)\ <C 1, which is the perturbativity condition 
previously given in Eq. (|B.6|) , 

At this stage, the result of the perturbative block diagonalization is: 

B \ B* M D M ) \ 0(S1 3 ) M'g j ' 

We can neglect the C(n 3 ) terms above. Thus, the remaining task is to re-diagonalize the two diagonal 
blocks above. However, as long as we work self-consistently up to second order in perturbation theory, 
no further re-diagonalization is necessary. Indeed, the off-diagonal elements of M'® and M'^ are of 
C(n 2 ). However, in the Takagi diagonalization, the off-diagonal terms of the diagonal blocks only effect 
the corresponding diagonal elements at 0(n 4 ) which we neglect in this analysis. The diagonal elements 
of .M^P and M'm also contain terms of 0(n 2 ), which generate second-order shifts of the diagonal 
elements relative to the my obtained at step [1]. These are easily obtained from the diagonal matrix 
elements of Eqs. ()B.13p and (|B.14p after making use of Eq. (|B.15P : 



mil 



, [ReByf] 2 [ImBiyf 2imf Re B Vj < lmBy r 

Tny + > < — — h — — 1 zzo =o } > l^- 1 '; 

^-^ I mn — mn mn 4- mn m , — m , 

(B.18) 



■ mji — my my + my m , — m -, 

j'=N'+l' (. 3 J 1 3 



2imy Re/3jy Im Bi'ji 




2 2 

m,-, - m% 
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with i' = 1', ..,7V' and j' = N' + l',..,N' + M'. Although the my are real and non-negative by 
construction, we see that the shifted mass parameters my are in general complex. Thus, to complete 
the perturbative Takagi diagonalization, we perform one final step. 

[3] The diagonal neutralino mass matrix is given by: 

M D N+M = ( M J * ) A7^t = diag(m?, • • • , < +M ,) , (B.19) 

where V is a suitably chosen diagonal matrix of phases 

V = diag(e-^i' , . . . , e -i **'+J*' ) , (B.20) 

such that the elements of the diagonal mass matrix A4® +M (denoted by m^) are real and non-negative. 
We identify the m V y with the physical neutralino masses. The unitary neutralino mixing matrix is then 
identified as: 

£). (B.21) 

Starting from Eqs. (fB~T7j) and (jB18|) . one can evaluate V to second order in the perturbation £1. In 
particular, for <C a, we have a + €\ + ie2 — (a + e\)e tt2 / a . From this result, we easily derive the 
second-order expressions for the physical neutralino masses m?/ : 

mf V L- ^ J + 1——LA-L L j' = i' . . . jv' b.22 

f ([Re^] 2 _[Img^] 2 1 [/ = iY ' + 1 ',..., iV ' + M']. (B.23) 

and the phases 



m -i ~ m,/ 



., ^ tn%> — Triji m^ + my 



N'+M' 

<k' - - J2 - <-2 -2 \ Rei % Im %' . [*' = !'> ■ ■ • . N '] > (B.24) 

N' _ 

- 52 _ -ax Re %' Im %' > [/ = iV' + 1' , . . . , iV' + M'] , (B.25) 

^ rnj>(rni,-m f ) 

This completes the perturbative Takagi diagonalization of the mass matrix for JV-dimensional and M- 
dimensional subsystems of Majorana fermions weakly coupled by an off-diagonal perturbation. As noted 
in Eq. (|B.6p . the perturbation theory breaks down if any mass from the A^-dimensional subsystem 
is nearly degenerate with a corresponding mass my from the M-dimensional subsystem (assuming that 
the corresponding residue, Re23jy, does not vanish). We provide an analytic approach to this case of 
near-degeneracy in Appendix C. 
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Appendix C: Degenerate mass eigenvalues 

If the value of one of the diagonal Ai^ elements, my, is nearly equal to one of the diagonal A4 M 
elements, say rn£>, and the corresponding residue He Bye does not vanish [cf. Eqs. (|B.22|) and (|B.23[) ] . 
then the techniques for degenerate states must be applied to diagonalize the full (N + M) x (N + M) 
matrix. We begin with the matrix Mn + m given in Eq. (|B.5p . which contains off-diagonal blocks of 
O(X), which characterizes the small couplings between the original MSSM matrix and the new USSM 
singlino/gaugino submatrix. 

We first interchange the first row and the k'th row of A^tv+m followed by an interchange of the 
first column and the k'th column, in order that my occupy the element of the matrix^ Next, we 
interchange the second row and the ^'th row followed by an interchange of the second column and the 
£'th. column, in order that occupy the 2"2" element of the matrix. This sequence of interchanges 
has the effect of grouping the two nearly degenerate diagonal elements next to each other, resulting in 
a new matrix M.' N+M with the following structure: 

\ 



M 



N+M 







( rriin 5 


A 


5 mi" 






■M-N+M-2 



(C.l) 



where the parameter 5 and the submatrix A are of O(X). The submatrix A4n+m-2 is no longer 
diagonal, although its new off-diagonal elements are all of 0{X). Thus, we may perform a perturbative 
Takagi diagonalization using the block diagonal unitary matrix, diag(l2x2 , jy N + M ~ 2 ^ with 

M N+ M-2 = (N N+M ~ 2 y M n +m-2 (N N+M -y = dmg(m' 3 „,m' A „ . . . ,m N „ +M „) , (C.2) 

where the m'-„ [j" = 3", 4", . . . , N" + M"\ are slightly shifted from the original non-degenerate {mj»}, 
{nTjii} by the perturbation of 

As a result of this procedure, the matrix M N+M in Eq. (|C.ip is modified by replacing the submatrix 
■M-N+M-2 by a diagonal matrix with perturbatively shifted diagonal elements, M.n+m-2- The off- 
diagonal blocks A and A T , are perturbatively shifted as well, but these shifts can be neglected as these 
effects are of higher order in the perturbation X. We denote the resulting matrix by .M jy+M • 

The complex parameter 5 couples the two near-degenerate states with mass parameters mi» and 
TO2». By definition of near-degeneracy, \m\n —myl <C S, so one cannot use perturbation theory in 
5 ~ 0{X). Instead, we shall perform an exact Takagi diagonalization of the 2x2 block ( m j" ^ f/ ^ of 



22 To distinguish the ordering of the physical neutralino states that arises from the manipulations performed in this 
appendix from the ordering of states based on Eqs. (|B.3[) and (|B.4[) . we employ double-primed subscripts here. 

23 For consistency with the second-order perturbative results of Appendix B, this diagonalization should be carried out 
including all contributions quadratic in X. 
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M-n+Mi using the results of Appendix A. 2: 



W* \ „ / 

oT 1 M N+M oT i 



m' i; , 
m' 2 n 




A T W^ 


■M-N+M-2 



(C.3) 



where the elements of the 2x2 unitary matrix W (which is denoted by W in Appendix A) can be 
determined in terms of Tnyi , and 5 using the formulae of Appendix A. 2. The (non-negative) diagonal 
masses vn-^u and m' 2 H are obtained from Eq. (|A.18|) : 

i' v , 2 „ = jmf„ + m|/ + 2\5\ 2 =f yj (m 2 ,,, - m 2 ,,) 2 + 4\m v ,5 + m 2 »5*| 2 j . (C.4) 



mi 



Note that if S is real, the quantity under the square root is a perfect square, in which case Eq. (|C.4|) 
reduces to the well-known expression: 

1 



2 L 



miii 



+ ni2" =F V (my - min) 2 + 45 2 



for real 5 . 



(C.5) 



If 5 is very small, the trajectories of the two eigenvalues nearly touch each other when the parameter 
M[ moves through the cross-over zone. A non-zero 5 value prevents the trajectories from crossing, 



keeping them at a distance > 5. In the 4'-5' zone, 5 is of first order in the ratio u/Msusy- in contrast, 
in the 2'-6' zone, 5 vanishes at first order due to the fact that V^g, ~ V^ 2 , ~ 0. However, as discussed 
below Eq. (|3.39p . these matrix elements acquire small non-zero corrections at higher order in v/Msvsy- 
Thus, we have two very different behaviors for 5, leading to the characteristically different evolution of 
the trajectories. These two cases are illustrated by the dashed lines in the two panels of Fig.O on the 
left for 5 — > and on the right for moderately non-zero 5 values. 

We may now apply the perturbative block diagonalization technique of Appendix B to complete the 
Takagi diagonalization of Eq. (|C.3p . The effect of this step is to shift the diagonal masses at second 
order as indicated in Eqs. (|B.22|) and (|B.23p . We finally arrive at the physical neutralino masses: 

N"+M" 



ph — / . 

m\„ ~ m v , + 



E 



[Re(W*A)i» f 



] 2 [Jm(W*A) 



m ■ 



m'-n + m'.„ 

1 J 



1" , 2"1 



(C.6) 



A' 



ph / 

m-„ ~ rrijii 



I [Re(W*A), v f 

^— Tn 1 - - — Try... 



m- 



' Im (W* A) in f 
m! V i + m'-,, 



\j" = 3",4", ... , N" + M"\ . (C.7) 



Since the appearance of rfi'y, and m' 2 „ [given by Eq. (|C4[) ] takes care of the near-degeneracy via an exact 
diagonalization (within the near-degenerate subspace), the results for the physical masses given above 
provide a reliable analytic description. 



The sizes of the second-order perturbative shifts in Eqs. (|C.6p and (|C.7p vary with the parameter 
M[ as the m'-„ \j" = 3", .., (N" + M")\ depend on M{. The effect of these shifts can be discerned in the 
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Figure 9: The evolution of the neutralino masses near the cross-over zone 1 (left) and near the crossover 
zone 2 (right) when varying the U(l)x gaugino mass parameter M[. The red dashed lines represent the 
masses of the diagonalized 2x2 matrix and the black solid lines after the subsequent approximate 
diagonalization of the full 6x6 matrix [Eq. IC.3|) A 



two cases considered above — in the cross-over zone 2'-6' with very small S, and in the cross-over zone 
4'-5' with moderately small S, as shown by the solid line trajectories of Fig.[9j 

Thus, we have demonstrated that an analytic perturbative treatment of the neutralino mass matrix 
can be carried out, and all of its features understood, even in the case of a pair of near-degenerate states. 
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